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[Text]  Abstract 

Research  on  an  RF  superconducting  cavity  at  Beijing  began  in  1988.  After  three 
years  of  hard  work,  substantial  progress  has  been  obtained,  and  has  been  well 
recognized  by  experts  worldwide.  The  preparation  and  procedures  for  the  low- 
temperature  superconducting  experiment  are  described  in  detail.  The  paper  also 
summarizes  progress  and  accomplishments  in  RF  (radio-frequency)  superconduct¬ 
ing  experimental  technique ,  microwave  and  phase— locked  measurement ,  post- 
purification  of  cavity,  computer  simulation  design  and  computer  control,  and 
data  acquisition  and  processing.  Presently,  an  accelerating  gradient  of  8.6 
MV/m  (Qg  “  6.5  X  10®)  has  been  obtained  in  a  1.5  GHz  niobium  cavity  operating 
in  the  CW  mode.  At  2  K,  it  was  measured  that  Qg  =  8  x  10®.  These  experimental 
results  are  at  the  cutting  edge  in  the  international  field. 

I .  Introduction 

RF  superconducting  technology  has  been  used  in  accelerators  for  some  time.  In 
1965,  a  lead-plated  superconducting  cavity  was  developed  at  Stanford 
University  to  accelerate  electrons.  After  more  than  two  decades  of  effort, 
substantial  progress  has  been  made  to  improve  the  shape,  material  performance 
and  accelerating  gradient  of  superconducting  cavities.^  Presently,  more  than 
20  laboratories  worldwide  are  using  or  building  high-energy  superconducting 
accelerators  for  heavy  ions  and  electrons. 

The  characteristics  of  a  superconducting  cavity  include  high  gradient,  large 
beam  diameter,  capability  of  operating  in  the  CW  mode,  and  low  operating 
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frequency.^  There  are  advantages  iti  using  a  superconducting  cavity  in  a  high- 
energy  intense  electron  beam  (e-beam)  accelerator.  A  superconducting 
accelerator  is  particularly  appropriate  for  a  free  electron  laser  (FEL) 
because  it  has  very  stringent  requirements  for  beam  quality,  intensity,  and 
mean  power.  In  the  late  80' s,  a  scheme  for  using  a  superconducting  cavity  and 
an  optical  cathode  high- intensity  injector  was  presented  with  the  objective  of 
increasing  the  brightness  of  an  FEL  or  TEV  linear  collider.^ 

Today's  superconducting  cavity  has  excellent  RF  characteristics.  In  the 
laboratory,  a  single  cavity  has  an  accelerating  gradient  of  25  MV/m  (Qq  «  10®). 
Based  on  BCS  theory,  the  accelerating  gradient  of  a  pure  niobium  (Nb)  cavity 
can  be  as  high  as  50  MV/m.^  The  experimental  value  is  still  far  less  than  that 
value.  Through  a  large  number  of  experiments  and  continuous  investigation,  it 
is  recognized  that  the  surface  characteristics  of  the  superconducting  cavity 
are  critical  to  its  overall  performance.  Processes  such  as  multiple  e-beam 
smelting,  increasing  RRR  (residual  resistance  ratio),  chemical  washing,  high- 
temperature  treatment  in  ultrahigh  vacuum,  and  artificial  oxide-film  formation 
can  effectively  improve  its  performance.  Nevertheless,  field-effect  emission 
is  still  the  limiting  factor  preventing  the  acceleration  gradient  from  rising 
further.  Some  large  laboratories  (e.g.,  KEK,  DESY,  CERN,  and  Cornell)  are 
looking  for  new  methods  and  techniques  to  modify  the  cavity  surface  in  order 
to  find  a  breakthrough  in  acceleration  gradient . 

The  heavy-ion  superconducting 
cavity  group  at  Beijing  Univer¬ 
sity  has  done  a  great  deal  of 
work  on  post-treatment  of  the 
cavity  to  raise  the  acceleration 
gradient.  As  a  result,  a  series 
of  effective  measures  has  been 
established.  Substantial  prog¬ 
ress  has  been  made  in  a  rela¬ 
tively  short  period  of  time.  We 
have  measured  a  gradient  of  8.6 
MV/m  (in  CW  mode).  At  2  K,  a 
cavity  quality  factor  of  8  x  10® 
has  been  obtained.  This  paper 
presents  methods  and  techniques 
used  to  improve  [the  quality 
factor]  Q  and  the  acceleration  Figure  1.  Elliptical  Cavity 

gradient 

II.  Design  Analysis  and  PC -Based  Program 
1.  Design  Analysis 

When  an  Nb  cavity  is  running  in  a  superconducting  state,  its  resistance  (R^) 
is  of  the  order  of  nO.  It  is  approximately  5  orders  of  magnitude  lower  than 
the  Rg  of  a  copper  cavity  at  liquid-nitrogen  temperature.^  This  extremely  low 
resistance  is  a  unique  feature  for  Nb.  An  important  point  is  that  the  shape  of 
the  superconducting  cavity  has  very  little  impact  on  the  RF  power  conversion 
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efficiency.  Hence,  cavity  design  should  address  other  issues  such  as  effects 
due  to  dark  arc  and  higher-order  modes  (ROMs).  Figure  1  shows  an  elliptical 
superconducting  cavity.  The  beam  tube  has  coupling  openings  and  both  ends  have 
cutoff  waveguides.  A  portion  of  the  accelerating  gradient  has  a  transitional 
circular  arc  and  it  also  has  a  large-aperture  beam  tube  to  effectively 
suppress  the  dark-arc  effect  and  beam-induced  field.  The  base  mode,  HOMs,  and 
wake— field  effects  of  a  1.5  GHz  LE— type  elliptical  single  harmonic  cavity 
(^=1)  have  been  analyzed  using  programs  such  as  SUPERFISH,  URMELT,  and  TBCI. 
Relevant  results  have  been  published  elsewhere. Figure  2  shows  the  wake 
field  caused  by  a  single  focused  pulse.  An  HP  8753C  network  analyzer  was  used 
to  measure  all  the  modes  below  3.0  GHz  associated  with  a  1.5  GHz  Nb  cavity  at 
room  temperature;  the  results  are  essentially  consistent  with  the  computer 
analysis.  Table  1  is  a  comparison  of  measured  and  calculated  values. 


Figure  2.  Wake  Field  Simulation  of  Elliptical  Cavity 


Table  1 .  Comparison  Between  Measurement  and  Computation 


MODE 

URMELT 

Me  as . 

MODE 

URMELT 

Me  as . 

(MHz) 

(MHz) 

(MHz) 

(MHz) 

TMOIO 

1470.22 

1467.70 

TMllO 

2058.45 

2043 . 54 

TMOll 

2792.13 

2800.40 

TMlll 

2621.95 

2585.17 

TM020 

2977.62 

2975.84 

TM112 

2955.58 

2974.58 

TElll 

1889.80 

1868.32 
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2.  Running  Software  Programs  on  PC 

In  computer  design  analysis,  an  outstanding  accomplishment  was  achieved. 
Programs  such  as  URMELT  and  TBCI  have  been  successfully  transplanted  onto  a 
386-PC.  To  run  a  large  program  on  a  PC,  the  first  issue  is  to  solve  the 
internal  memory  problem.  The  code  and  data  segment  of  programs  such  as  URMELT 
and  TBCI  are  much  larger  than  640  K,  which  makes  it  impossible  to  run  directly 
in  a  DOS  environment.  In  the  past,  the  conventional  method  has  been  to  use 
overlays  for  the  code  and  to  open  temporary  files  for  the  large  data  segment. 
However,  this  requires  the  careful  analysis  of  each  module  and  every 
executable  statement  in  the  source  code.  It  takes  an  incredible  amount  of 
work. 

Our  approach  is  to  use  virtual  memory;  this  divides  up  a  portion  of  the  hard 
drive  for  virtual  memory.  When  a  386-PC  operates  under  a  protected  virtual 
address  mode,  it  is  capable  of  managing  16  M  of  real  memory  and  64  T  (T  =  M  x 
M)  of  virtual  memory.  This  makes  it  possible  to  run  large  programs  on  a  PC.  In 
this  case,  the  size  of  the  program  is  determined  by  the  capacity  of  the  hard 
drive.  On  our  HP-312  PC,  with  10  M  internal  memory  (2  M  real  and  8  M  virtual), 
both  URMELT  and  TBCI  could  be  run  normally.  The  speed  is  approximately  one- 
third  to  one-half  of  that  on  the  VAX.  However,  with  multiuser  time-sharing  on 
a  VAX,  it  is  comparable  in  speed.  Figure  2  shows  some  TBCI  results  obtained 
from  the  HP-312  PC. 

The  fact  that  these  programs  can  run  on  a  386-PC  is  extremely  convenient  for 
the  research  group,  not  only  saving  us  time  but  also  money. 

III.  Cryostat  and  Ultrahigh-Vacutam  System 

1 ,  Cryostat 

The  1.5  GHz  Nb  superconducting  cavity  experiment  has  to  be  performed  at  2-4.2 
K.  To  this  end,  a  cryostat  was  constructed.  After  two  years  of  effort,  we 
completed  the  fabrication  and  installation  of  a  system  consisting  of  a  large 
cryostat,  a  100-liter  liquid  helium  (He)  Dewar,  two  100-liter  liquid-nitrogen 
(N2)  Dewars,  a  liquid-He  duct,  an  He  recovery  duct  and  gauges  for  temperature 
and  liquid-He  height.  The  construction  of  this  apparatus  ensures  the  smooth 
execution  of  superconductivity  experiments . 

The  100-liter  cryostat  is  the  major  equipment  required  to  perform  the 
experiment.  The  development  of  such  a  large  low-temperature  liquid-He  device 
had  never  been  reported  domestically.  Due  to  lack  of  information,  the  initial 
design  was  sent  to  DESY  (Germany)  for  review  by  foreign  experts.  Our  German 
colleagues  came  back  with  10  major  modifications.  The  design  was  then  revised 
and  made  simpler  and  more  practical.  The  new  design  stresses  safety  and  the 
ultrahigh-vacuum  requirement  associated  with  the  Nb  superconducting  cavity. 
Figure  3  is  a  sketch  of  the  cryostat,  which  is  approximately  175  cm  high  with 
a  maximum  outer  diameter  of  75  cm.  The  Nb  cavity  and  the  cryostat  are  on  two 
separate  vacuum  systems.  The  cryostat  uses  a  high-speed  molecular  pump  and  the 
system  is  required  to  provide  10“^-10"^  Pa  of  vacuum.  The  Nb  superconducting 
cavity  has  an  ultrahigh-vacuum  system. 
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After  about  a  year  of  operation,  the 
entire  cryostat  system  has  essentially 
met  the  design  requirements.  Liquid  He 
in  the  cryostat  evaporates  at  a  rate 
of  1  percent  per  hour.  The  100-liter 
liquid-He  Dewar  evaporates  at  a  rate 
of  less  than  1.5  percent  per  day. 


2.  Ultrahlgh-Vacuum  System 

The  internal  surface  of  the  supercon¬ 
ducting  cavity  must  be  extremely 
clean.  Otherwise,  it  has  a  significant 
impact  on  Q  and  Eacc-  The  1-5  GHz  Nb 
cavity  is  vertically  mounted  in  the 
cryostat  and  is  connected  to  the 
vacuum  line  in  an  overhead  manner  to 
prevent  any  debris  from  falling  into 
the  cavity.  Both  ends  of  the  cavity 
are  sealed  with  an  Indium  (In)  sealing 
surface.  The  inner  diameter  of  the  In 
seal  flange  must  be  designed  according 
to  the  actual  size  of  the  beam  pipe 
flange.  It  must  be  highly  accurate. 

The  diameter  of  the  sealing  In  wire  is 

1  mm.  These  two  large— aperture  sealed 
connections  can  directly  influence  the 
ultrahigh-vacuum  achievable  by  the 
system.  Especially  when  liquid  He  is 
in  a  superfluid  state,  these  connections  are  the  weak  links.  Two  cryostatic 
absorption  pumps  and  an  ion  pump  provide  an  oil-free  pumping  system  for  the  Nb 
cavity .  The  entire  system  is  metal— sealed  with  ultrahigh-vacuum  metal  valves . 
The  RF  feed-throughs  are  50-0  ceramic— seal  coaxial  connectors  (total  99).  An 
ultrahigh-vacuxim  jewel  valve  is  installed  at  the  end  of  the  vacuum  line  as  a 
control  valve  for  helium  treatment  experiments. 

This  system  can  reach  10“®~10~^Pa  of  vacuum  before  using  liquid  N2  to  lower  the 
temperature.  After  adding  liquid  N2  and  He,  it  can  reach  10“®Pa.  After  repeated 
tests  with  liquid  He,  the  system  was  found  to  be  reliable  for  low-temperature 
physical  experiments.  The  first  phase  of  work  was  thus  completed  smoothly, 

3.  Low-Pressure,  Low-Temperature  Experiment 

The  Q  value  of  a  superconducting  cavity  varies  with  temperature.^  At  4.2  K,  Q 
is  (2-4)  X  10®  for  the  1.5  GHz  Nb  cavity.  At  2  K,  Qq  -  (6-9)  x  10®.  Hence,  it 
is  essential  to  achieve  an  ultralow  temperature  in  order  to  reach  a  higher 
gradient.  In  order  to  reach  an  ultralow  temperature,  we  employed  a  method  to 
lower  temperature  and  reduce  pressure  by  a  mechanical  pvimp.  After  many  tests, 

2  K  was  finally  achieved  in  the  100-liter  cryostat  and  Qq  measured  8  x  10®.  At 
2.12  K,  liquid  He  is  a  superfluid,  posing  a  serious  challenge  to  the  vacuum 
components  in  the  liquid-He  tank  because  superfluid  He  is  extremely  permeable. 
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After  niimerous  low-temperature,  low-pressure  experiments,  the  vacuum  inside 
the  Nb  cavity  remained  unchanged.  This  demonstrates  the  integrity  of  the  In 
seals  and  other  metal  seals.  It  is  rare  in  China  for  such  a  large  liquid-He 
vessel  to  reach  such  a  low  temperature . 

4.  Shielding  of  the  Earth's  Magnetic  Field 

The  Earth's  magnetic  field  has  a  great  deal  of  impact  on  the  RF  characteris¬ 
tics  of  the  superconducting  cavity.  Hence,  it  is  necessary  to  establish  a 
magnetic  shield  around  the  Nb  cavity.  We  added  several  layers  of  thin  film  of 
a  high-/i  material,  which  is  an  amorphous  metal  (Fe-Ni-Si-B) ;  the  thin  ribbon 
is  0.04  mm  x  100  mm  (thickness  x  width),  total  weight  6  kg.  The  axial  and 
horizontal  magnetic  field  inside  the  cryostat  were  measured  by  technicians 
from  the  Department  of  Geophysics  at  Beijing  University  using  a  model  DM  2220 
Digital  Magnetometer  (U.S.-made).  It  was  found  that  the  axial  magnetic  field 
is  one-seventh  and  the  horizontal  magnetic  field  is  one-twentieth  of  the 
Earth's  magnetic  field.  This  essentially  satisfies  the  requirement  for  the 
superconducting  experiment . 

IV.  Post-Treatment  of  Nb  Cavity 

After  a  superconducting  Nb  cavity  is  machined  and  welded,  it  needs  to  be 
polished  mechanically,  electrically,  and  chemically  and  then  purified  after 
putting  on  a  titanium  (Ti)  film.  As  a  result  of  these  processes,  a  high- 
quality,  high-thermal-conductivity  inner  surface  is  produced.  These  processes 
help  improve  the  quality  factor  Q  and  accelerating  gradient  Egcc-  This  involves 
materials  science,  chemistry  and  thermal  treatment.  In  order  to  obtain  the 
ideal  RF  superconducting  conditions,  the  following  work  was  done. 

1.  Analysis  of  Nb  Material 

The  superconducting  cavity  inside  an  electron  accelerator  is  usually  made  of 
high-purity  thin  Nb  plate.  Commercial-grade  Nb  cannot  be  used  directly; 
several  times  of  e— beam  smelting  are  required  to  raise  its  purity.  After 
raising  purity,  it  is  necessary  to  improve  its  thermal  conductivity.  The 
normal  method  is  to  evaporate  a  Ti  coating  on  the  Nb  surface  in  order  to 
remove  the  impurity  in  Nb  by  way  of  the  combining  process. 

We  analyzed  and  compared  the  high-RRR  Nb  supplied  by  DESY  to  the  e-beam- 
smelted  Nb  processed  by  the  Institute  of  Nonferrous  Metals .  Table  2  shows  the 
semiquantitative  chemical  composition  of  the  two  specimens.  Sample  #1  is  from 
DESY  and  is  not  very  impressive  from  the  standpoint  of  impurity.  Figure  4,  an 
electron  spectrum  of  the  specimen,  shows  that  the  surface  is  covered  by  a 
layer  of  lhj205  several  hundred  angstroms  in  thickness.  After  repeated  RF 
superconducting  experiments,  it  was  verified  that  this  oxide  layer  could 
stabilize  the  RF  characteristics  of  the  cavity.  Presently,  this  high-quality 
oxide  film  is  being  added  to  the  purification  process  at  KEK  and  CEBAF  as 
well. 
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Table  2.  Semi-Quantitative  Chemical  Analysis  of  Two  Nb  Specimens 
Sample  #1  [Figures  in  percent] 


Ni 

Zr 

Rh 

Sn 

T1 

Ta 

Mg 

Pt 

Ti 

0.001 

0.003- 

0.01 

0.01 

0.001 

0.001 

0.1 

0.001 

0.002 

Sample  #2  [Figures  in  percent] 


Hf 

Zr 

Tl 

A1 

Pt 

Cu 

Ti 

0.01 

0.003 

0.001 

0.003 

0.001 

0.001 

0.003 

Sample  *1 


X00542  8-3-89 

REGIOM  I  IHMi 

KPS  RL  CREjSO  S.NIS7  96  SEC 


0.050  STEP 


2.  Chemical  Polishing 

The  1.5  GHz  Nb  cavity  manu¬ 
factured  by  Dornier  Corpo¬ 
ration  is  made  of  high-pu- 
rity  Nb  with  RRR  *  300. 

Usually,  the  thermal  break¬ 
down  threshold  is  higher 
when  RRR  is  high.  It  will 
be  easier  to  attain  a  high¬ 
er  accelerating  gradient 
EaccI  however,  it  is  more 
difficult  to  chemically 
clean  this  high-RRR  Nb.  In 
the  acid  washing  process, 
it  is  easy  to  incorporate 
H'*’  or  phosphate  into  the  Nb 
surface  to  cause  substan¬ 
tial  deterioration  of  sur¬ 
face  characteristics.  In 
order  to  avoid  this  situa¬ 
tion,  we  conducted  numerous  chemical  polishing  experiments  and  finally 
established  an  effective  polishing  method. 


B  »  E  (cV) 

Figure  4.  Photoelectron  Spectrum  of  Nb 
Surface 


Chemical  polishing  was  done  according  to  the  method  used  at  DESY  and  KEK.^^ 
Figure  5  [photos  not  reproduced]  shows  scanning  electron  micrographs  of  a 
specimen  before  and  after  chemical  polishing.  It  is  obvious  from  the  pictures 
that  the  crystals  are  well  aligned  on  the  Nb  surface  after  chemical  polishing. 
The  polishing  solution  is  composed  of  a  fixed  ratio  of  HF  (40%),  HNO3  (65%), 
and  H3P0^  (85%)  .  The  number  of  washes  and  the  time  of  exposure  must  be 
carefully  controlled.  The  Nb  cavity  needs  to  be  cooled  during  the  process. 
Cleaning  after  acid  polishing  is  very  critical.  The  resistivity  of  the 
deionized  water  used  must  be  monitored  frequently;  especially  in  the  final 
rinse,  it  should  be  at  approximately  18  MQ.  The  cavity  is  mounted  on  a  special 
support  during  acid  polishing.  It  can  be  rotated  to  allow  the  acid  to  come 
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into  contact  with  the  cavity  wall  evenly  and  to  make  the  discharge  of  the 
solution  easier.  Acid-resistant  flanges  and  valves  are  used  on  both  sides  of 
the  beam  pipe.  The  first  chemical  polishing  prior  to  thermal  treatment  must 
remove  a  surface  layer  20-30  /i  in  thickness.  After  thermal  treatment,  it  needs 
to  be  chemically  polished  one  more  time  to  remove  another  5-10  /i. 

3.  High-Temperature  Treatment 

We  began  conducting  RF  superconducting  experiments  in  October  1990.  In  the 
first  two  experiments,  only  chemical  polishing  was  done  and  the  results  were 
not  satisfactory.  After  reviewing  the  data  obtained  in  KEK,  DESY,  and  Cornell, 
our  preliminary  conclusion  was  that  stress  in  the  cavity  had  to  be  relieved. 
This  stress  comes  from  machining  and  welding.  According  to  a  1989  Cornell 
publication,  it  was  reported  that  the  Qq  for  their  1,5  GHz  Nb  cavity  was  two 
orders  of  magnitude  lower  than  the  expected  value  and  could  not  be  improved 
until  the  stress  was  relieved  by  heat  treatment.  DESY  also  presented  the  same 
issue  in  a  1989  meeting.  In  reference  to  their  experience,  a  thermal  treatment 
procedure  was  formulated  as  follows:  The  heat-resistant  is  done  in  a  vacuum 
furnace  with  a  vacuum  better  than  1  x  10“^Pa  and  at  a  temperature  of  850-900°C 
for  one  hour.  The  vacuum  pumps  should  not  be  turned  off  during  cooling.  The  Nb 
cavity  is  placed  in  a  Ti  box  during  the  heat  treatment.  Impurities  in  Nb  are 
evaporated  from  the  surface  at  high  temperature  and  may  be  readsorbed  during 
cooling.  With  a  Ti  box,  impurity  ions  are  easily  adsorbed  by  Ti  to  keep  the  Nb 
cavity  from  secondary  contamination. 

The  Nb  cavity  used  in  our  third  liquid-He-temperature  experiment  was  thermally 
treated;  this  is  the  very  cavity  with  which  we  had  a  major  breakthrough.  The 
measured  Qq  was  of  the  order  of  10®  (at  4.2  K)  .  With  reduced  pressure  and  lower 
temperature,  Qq  was  found  to  be  of  the  order  of  10®.  These  numbers  are 
consistent  with  results  reported  in  other  countries  and  they  are  also  in 
agreement  with  theoretical  values.  Furthermore,  we  also  noticed  that  the  heat- 
treated  cavity  was  very  reproducible.  We  did  not  see  Qq  drop  due  to  temperature 
excursion. 

4.  Ultraclean  Installation 

In  order  to  prevent  contamination,  the  cleaned  cavity  must  be  installed  in  an 
ultraclean  environment.  After  chemical  polishing,  the  cavity  is  filled  with 
ultrapure  water  to  prevent  any  impurity  from  entering  the  cavity.  It  is  then 
immediately  sent  to  an  ultraclean  room  for  installation.  Prior  to  installa¬ 
tion,  the  cavity  is  rinsed  again  with  ultrapure  water.  The  ultraclean  room  is 
equipped  with  an  air  curtain  and  air  filters.  It  also  has  a  Class— 100  work¬ 
bench  for  two  to  ensure  that  there  are  fewer  than  100  dust  particles  more  than 
0.5  fim  in  size  per  cubic  foot  over  the  workbench.  The  ntimber  of  such  dust 
particles  in  a  cubic  foot  of  normal  air  is  210,000.  When  installing  the  main 
coupler,  pure  nitrogen  evaporated  from  the  liquid-nitrogen  Dewar  is  fed  into 
the  cavity  through  the  heating  tube  to  create  a  positive  pressure  inside  the 
cavity.  This  is  to  prevent  air  from  coming  in  contact  with  the  cavity  surface 
to  minimize  the  chance  of  any  impurities  entering  the  cavity.  After  the  cavity 
is  connected  to  the  main  coupler,  nitrogen  is  injected  from  the  main  coupler. 
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These  two  components  are  then  connected  to  the  vacuum  system.  During  this 
process,  a  positive  pressure  is  always  maintained  inside  the  cavity. 

After  turning  on  the  ion  pump,  the  entire  system  is  baked  for  approximately  20 
hours,  and  will  reach  the  desired  vacuum  in  a  week  or  so. 

V.  Microwave  and  Phase -Locked  Measurement  Technique 

Microwave  measurement  is  a  key  step  in  RF  superconducting  experiments.  Under 
superconducting  conditions,  the  resistance  Rg  of  the  1.5  GHz  Nb  cavity  is  less 
than  a  few  hundred  ohms,  which  is  five  orders  of  magnitude  lower  than  the  Rg 
of  a  copper  cavity.  Hence,  the  quality  factor  of  the  superconducting  cavity  is 
as  high  as  10®--10®.  The  half-power  bandwidth  of  its  resonant  curve  is  only  a 
few  Hertz  wide  and  the  field-decay  time  constant  is  30-300  ms.  These  proper¬ 
ties  make  the  measurement  of  RF  parameters  and  search  for  the  resonant  point 
very  difficult.  From  the  very  beginning,  we  recognized  that  microwave  measure¬ 
ments  involving  such  a  high-Q  superconducting  cavity  were  never  attempted  in 
China  and  it  would  require  the  appropriate  test  methods  and  instrumentation 
with  high  accuracy  and  stability. 

1 .  Instrument  Setup 


HP5470B 

Figure  6 .  Microwave  Phase-Locked  Loop 

6  shows  the  microwave  measurement  circuit  used  in  the  low— temperature 
superconducting  experiment,  and  essentially  shows  our  instrument  setup  as 
well,  'pie  synthesized  signal  generator  is  an  HP  8663A,  which  has  a  frequency 
stability  of  10  ^°/day.  The  output  signal  has  a  resolution  of  0.4  Hz  at  near 
1.5  GHz,  which  is  less  than  the  half— power  bandwidth  of  the  superconducting 
cavity.  The  signal  generator  is  capable  of  scanning  and  phase-locked  tuning. 


When  searching  for  the  superconducting  point  by  frequency  scanning,  a  signal 
generator  must  have  the  kind  of  high  resolution  and  stability  as  those 
described  above. 

The  bandwidth  of  the  four-channel  digital  oscilloscope,  an  HP  54503A,  is  DC- 
500  MHz.  The  horizontal  time  axis  has  501  points.  (The  HP-IB  remote  data 
acquisition  system  has  1,024  points).  The  vertical  voltage  axis  is  eight-bit. 
Its  accuracy  is  1/2®  -  0.4  percent.  (The  HP-IB  remote  data  acquisition  system 
is  16-bit  with  an  accuracy  of  1.5  x  10'^.)  It  can  automatically  detect  the 
maximum  and  minimum  signal.  Furthermore,  it  is  capable  of  measuring  the 
voltage  at  a  specific  time  by  remote  control.  It  also  can  determine  the  time 
at  which  the  signal  is  at  a  specific  voltage.  The  oscilloscope  also  has  memory 
for  storage. 

The  power  meter  is  an  HP  437B  and  its  measurement  range  is:  100  kHz-50  GHz  in 
frequency,  -70  dBm  to  +44  dBm  in  power,  and  0.1  percent  in  accuracy.  The 
frequency  counter  is  an  HP  5350B.  It  has  an  accuracy  of  1  Hz  near  1.5  GHz  and 
the  gate  is  200  ms  wide. 

The  HP  PC-312  controller  has  two  CPUs.  It  not  only  can  function  as  a  386 
microcomputer  but  also  can  be  connected  to  the  signal  generator,  oscilloscope, 
power  meter,  frequency  counter,  and  laser  printer  through  an  IEEE-488.2 
interface  to  form  a  fully  automated  data  acquisition,  processing,  storage,  and 
print-out  system. 

The  20— W  RF  solid-state  amplifier  is  the  final— stage  power  amplifier  of  the 
measurement  system.  It  was  constructed  by  Institute  54  of  the  Ministry  of 
Machine-Building  and  Electronics  Industry  (MMEI).  The  frequency  range  of  the 
amplifier  is  1.4  GHz-1.6  GHz.  When  Qq  of  the  Nb  cavity  is  of  the  order  of  10®, 
the  maximum  power  output  meets  the  experimental  requirement. 

The  system  employs  directional  couplers  such  as  HP  778D  (20  dB)  ,  Narda  4242-20 
and  TT3-W  (10  dB) ,  and  an  HP  8470B  crystal  detector  and  a  DMM-6-1500  dual 
balance  mixer  (DBM) . 

2.  Phase-Locked  Loop  (FLL) 

The  half— power  bandwidth  at  the  resonant  point  of  the  superconducting  cavity 
is  only  a  few  Hz.  However,  the  frequency  drift  caused  by  external  mechanical 
vibration  is  approximately  20-30  Hz.  The  frequency  drift  during  cooling  is 
approximately  0.1  MHz/K.  To  this  end,  the  microwave  measurement  system 
requires  a  PLL  technique  to  allow  the  frequency  of  the  signal  generator  HP 
8663A  to  track  the  frequency  change  of  the  cavity  in  order  to  ensure  that  the 
superconducting  cavity  remains  in  a  resonant  state. 

Figure  6  shows  the  PLL  measurement  circuit.  The  microwave  signal  generated  by 
the  HP  8663a  (10  mW)  travels  through  the  20-W  power  amplifier,  isolator,  and 
variable  attenuator  to  reach  the  HP  778D  directional  coupler.  At  that  point, 
the  main  signal  is  directly  fed  into  the  main  coupler  of  the  superconducting 
cavity  from  the  exit  of  the  directional  coupler.  In  addition,  a  -20  dB  signal 
is  extracted  from  the  incident  wave  to  be  transmitted  to  one  port  (LO)  of  the 
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DBM  after  passing  through  a  mechanical  phase  shifter.  Another  microwave  signal 
is  extracted  from  the  superconducting  cavity  to  be  transmitted  to  the  RF  end 
of  the  DBM.  A  phase-error  signal  will  be  sent  after  comparing  the  two  signals. 
This  error  signal  goes  through  a  PI  amplifier  and  is  then  transmitted  to  the 
FM  port  of  the  HP  8663A  to  form  a  stable  feedback— control  circuit. 

3.  Main  Coupler  and  Pick-Up  Probe 

The  microwave  power  is  fed  into  the  superconducting  cavity  through  the  main 
coupler.  The  pick-up  probe  extracts  the  microwave  signal  from  the  cavity  for 
measurement.  As  far  as  power  input  and  pick-up  signal  are  concerned,  the 
external  quality  factors,  Qg^t,  are  defined  as  follows: 


(input  portion) ,  =  2nU/Wp  (pick-up  portion)  (1) 

where  U  is  the  energy  stored  in  the  cavity,  and  Wg  and  Wp  represent  the  energy 
loss  per  cycle  at  the  internal  resistance  of  the  signal  generator  and  that  in 
the  signal  pick-up  circuit,  respectively.  It  is  known^  that  the  quality  factor 
with  load,  Q^,  is  related  to  [the  unloaded  quality  factor]  Qq  as  follows: 


l/!?i  =  l/Oo  + 


(2) 


In  order  to  express  the  degree  of  coupling, 
let  us  introduce  the  coupling  factors  for  the 
input  circuit  and  pick-up  circuit  as  and 
^2- 


Pi  =  Pi  »  Oo/Oext,  (3) 

Substituting  and  into  equation  (2),  we 
have : 


Qo  =  Qi  (1  +  ySi  +  ^2)  (4) 

The  quality  factor  of  the  cavity  with  load, 
Qi,  can  be  measured  experimentally.  The 
unloaded  quality  factor,  Qq,  can  only  be 
determined  after  fix  and  ^2  nre  obtained. 
Figure  7  is  a  schematic  diagram  of  the  main 
coupler.  The  Qq  value  is  approximately  10®- 
10®.  In  order  to  feed  as  much  microwave  power 
as  possible  into  the  cavity,  it  is  required 
that  of  bi'®  main  coupler  be  comparable 

to  Qq.  Since  of  the  pick-up  probe  is  two 

orders  of  magnitude  higher  than  Qq,  the 
magnitude  of  the  signal  picked  up  by  the 
probe  has  little  effect  on  the  measurement 
neglected  in  calculating  Qg. 


of  Qq.  ^2  of  the  probe  can  be 
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The  external  Q  values  are  designed  to  be:  O^xt^  “  (1~3)  x  10^,  Qext^  =  1  x  10^^ 
and  the  length  of  the  main  coupler  antenna  and  that  of  the  probe  are 
determined  accordingly.  These  values  were  first  estimated  using  a  computer  and 
then  adjusted  experimentally;  finally,  they  were  verified  during  superconduc- 
tivity.  If  necessary,  the  devices  will  be  reinstalled;  this  is  a  delicate 
piece  of  experimental  work.  At  liquid-He  temperature,  the  value  of  Q  varies 
due  to  thermal  expansion  and  contraction.  The  pick-up  probe  is  made  of  a  50-0 
hard  coaxial  cable.  During  cooling,  the  insulating  layer  between  the  inner  and 
outer  conductor  contracts  several  dozen  times  more  than  the  metal.  The  result 
is  that  its  is  two  orders  of  magnitude  higher  than  the  designed  value  at 
liquid-He  temperature.  We  encountered  similar  problems  with  the  main  coupler. 


The  antenna  of  the  main  coupler  is  a  piece  of  50-0  coaxial  cable.  There  is  a 
vacuum  between  the  inner  and  outer  conductor.  One  end  is  connected  to  the 
power  feed  line  through  a  50-0  coaxial  ceramic  seal.  This  ceramic  connector 
can  only  operate  down  to  20  K,  The  inner  conductor  is  made  of  molybdenum.  The 
antenna  and  the  inner  conductor  are  threaded  together.  Initially,  the  antenna 
was  made  of  copper.  After  the  experiment,  it  was  found  that  the  antenna  almost 
came  loose  and  there  was  only  one  thread  left  connected.  This  is  because 
copper  has  a  higher  thermal  expansion  coefficient  than  molybdenum.  Through 
experimentation,  the  design  was  modified  and  materials  were  changed,  in  order 
to  make  experimental  values  obtained  at  room  temperature  consistent  with  those 
obtained  at  low  temperature. 


The  electric  field  distribution,  E^, 
in  the  cavity  was  calculated  using 
URMELT.  Based  on  the  relationship 
between  probe  voltage,  electric  field 
in  the  cavity  and  probe  length,  it  is 
possible  to  estimate  the  length  of  the 
probe.  Figure  8  shows  the  experimental 
method  to  measure  the  lengths  of  the 
antenna  and  the  probe.  An  auxiliary 
antenna  is  set  at  one  end  of  the  cavi¬ 
ty  and  its  length  is  adjusted  based  on 
the  reflective  scan  signal  to  make 
^  1.  The  output  power  of  the  signal 

generator  is  measured  with  a  power 
meter.  Then,  the  pick-up  power  of  the 
main  antenna  at  the  other  end  of  the 
cavity  is  measured.  The  length  of  the 
antenna  is  adjusted  to  obtain  a  suit¬ 
able  power  pick-up  to  output  ratio. 
Thus,  the  desired  external  Q  values  c; 


Figure  8.  Principle  of  Q^xt  Measurement 


be  obtained.  In  our  experiments, 


Qext^  =  1  x  10^,  Qext^  =  2  X  10^°.  The  latter  is  one  order  of  magnitude  lower  than 
the  designed  value  because  the  phase  comparator  requires  the  probe  to  provide 


a  signal  higher  than  5  dBm. 
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VI.  RF  Superconducting  Experiments  and  Results 

RF  superconducting  experiments  began  in  October  1990.  The  first  experiment  was 
to  test  the  entire  low-temperature  system.  The  low-temperature  and  vacuum 
system  was  found  to  perform  well;  it  was  safe  and  reliable  to  operate.  The 
system  is  capable  of  pumping  the  Nb  cavity  and  the  cryostat  jacket  to  10"®Pa 
and  10  '‘Pa,  respectively,  prior  to  cooling.  After  two  days  at  liquid-Nj 
temperature,  the  temperature  of  the  Nb  chamber  reaches  160  K.  During  a 
superconducting  experiment,  liquid  He  is  transported  from  the  Dewar  to  the 
cryostat  in  a  pipe  until  the  entire  Nb  cavity  is  submerged.  At  the  same  time, 
the  recovery  pipe  is  turned  on  to  receive  He  in  the  low-temperature  workshop. 
Afterward,  the  liquid-He  pipe  is  removed.  The  temperature  of  the  Nb  cavity  is 

4.2  K.  The  level  of  liquid  He  is  above  the  upper  flange  of  the  beam  pipe.  The 
temperature  and  liquid-He  level  are  monitored  with  an  Fe-Rh  thermistor  and  a 
liquid-He-level  meter,  respectively.  A  15-liter  mechanical  pump  is  attached  to 
the  cryostat  with  a  three-way  connector.  Its  exhaust  is  connected  to  the 
return.  The  reduced-pressure,  low-temperature  experiment  can  proceed  by 
turning  on  the  mechanical  pump  and  controlling  the  valve.  The  lowest  tempera¬ 
ture  achievable  is  2  K. 

1.  RF  Superconducting  Results 

The  RF  superconducting  experiment  is  designed  to 
investigate  the  RF  characteristics  of  the  1 . 5 
GHz  Nb  cavity  in  the  superconducting  state.  This 
kind  of  L— band  resonant  cavity  is  used  in  a 
linear  accelerator.  Experimentally,  the  unloaded 
Q  value  was  measured  as  a  function  of  accelera¬ 
tion  gradient,  Egcc.  and  temperature,  T.  Tempera¬ 
ture  was  varied  by  cooling  under  reduced  pres¬ 
sure.  Q  values  of  the  cavity  at  6.5  K,  at  4.2  K, 
and  at  several  points  between  4.2  K  and  2  K  were 
measured.  Such  measurements  were  made  over  and 
over  again  with  excellent  reproducibility. 

Figure  9  shows  a  curve  plotting  surface  resis¬ 
tance  Rg  against  T^/T.  The  relationship  between 
Rg  and  Qq  is  G  =  Rg  x  Qq,  where  G  is  the  geomet¬ 
ric  constant  of  the  cavity.  For  this  cavity, 

G  -  260  0.  Tg  is  the  critical  temperature  and  is 

9.2  K  for  Nb.  The  solid  line  in  Figure  9  is  the 
expected  Rg  versus  temperature  curve  derived  from 
BCS  theory.  The  triangles  are  experimental  points.  The  theoretical  values  are 
in  very  good  agreement  with  experimental  points  in  most  areas.  At  very  low 
temperatures,  the  experimental  points  deviate  from  the  theoretic  curve  and 
gradually  approach  a  fixed  value.  This  is  normal  and  the  value  is  known  as  the 
residual  resistance,  Rges,  of  the  superconducting  cavity.  This  is  also  the 
difference  between  RF  superconductivity  and  dc  superconductivity.  The 
magnitude  of  Rpes  is  dependent  on  Nb  surface  characteristics . ^ 

Figure  10  shows  the  Qq  versus  Eggc  curve.  A  power  experiment  was  performed  at 

4.2  K  and  2.25  K.  Qq  and  Eggc  were  measured  at  various  input  power.  As  shown  in 


TJT 

Figure  9.  Curve  of  Rg  Vs. 

Tc/T 
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the  figure,  Qg  decreases  with  in¬ 
creasing  Eacc-  The  rate  of  decline 
and  the  presence  of  an  obvious 
turning  point  are  criteria  to  judge 
the  RF  characteristics  of  a  cavity. 

The  Qo  versus  Egcc  curve  should  be 
relatively  flat  for  a  good  cavity. 

At  a  higher  gradient  of  acceler¬ 
ation,  Qo  should  drop  or  show  an 
apparent  turning  point.  No  turning 
point  was  found  in  our  experiment. 

However,  Qg  did  fall  somewhat.  In 
our  initial  power  experiment,  a 
gradient  of  acceleration  of  3  MV/m 
was  obtained.  After  repeatedly 
improving  techniques  such  as  chemical  polishing,  thermal  treatment,  reduced 
pressure  cooling,  and  microwave  PLL,  it  has  been  increased  to  8.6  MV/m.  This 
value  is  already  at  a  world-class  level.  Currently,  we  are  preparing  to 
conduct  the  more  difficult  He  experiment  to  further  improve  the  Qg  versus  Egcg 
curve  to  obtain  a  larger  Qg  at  higher  Egcc- 

In  the  figure,  Egcc  is  calculated  from  the  following: 

^acc  =  (\^/o.i),/io2”jr^"jr74p7Ti^^Fnr^  (s) 

where  0.1  is  the  effective  acceleration  length  of  the  cavity  (m) ,  (102  x  Qq) 
is  the  impedance  of  the  cavity  (Q) ,  Pf  is  the  incident  RF  power  (W)  and  Eacc 
in  units  of  (V/m)  .  Values  of  p  and  Qq  are  discussed  in  detail  in  the  following 
section. 

2.  Microwave  Measurement  Technique  Under  Superconducting  Conditions 

As  mentioned  earlier,  the  unloaded  quality  factor,  Qq,  of  the  1,5  GHz  Nb  cavity 
is  10®-10^  and  the  half-power  bandwidth  is  only  a  few  Hz.  The  field-decay  time 
constant,  however,  is  often  as  long  as  30-300  ms.  This  makes  the  measurement 
of  RF  parameters  and  location  of  the  resonant  point  very  difficult.  Further¬ 
more,  because  our  experimental  cavity  does  not  have  a  frequency  tuning 
mechanism,  it  makes  it  more  difficult  to  determine  the  resonant  frequency. 

(1)  Locating  the  Resonant  Point  by  Frequency  Scanning 

Because  the  HP  8663A  synthesized  signal  generator  has  a  very  high  frequency 
stability  and  a  0.4  Hz  frequency  scan  resolution,  precision  scanning  is 
carried  out  in  every  100  Hz  region  between  1.46  GHz  and  1.47  GHz  at  a  rate  of 
100  ms/Hz.  After  repeated  precision  frequency  scanning,  the  superconducting 
resonant  frequency  was  captured.  After  this  frequency  is  obtained,  the 
generator  is  switched  from  scanning  to  CW  mode.  In  addition,  the  PLL  begins  to 
operate.  Under  the  control  of  the  PLL,  the  superconducting  cavity  operates 
steadily  in  a  resonant  state. 


Figure  10.  Qq  Vs.  Eacc 
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(2)  Measurement  of  ^  by  Reflection 

The  value  of  ^  is  measured  by  using  the  transient  response  of  the  resonance 
circuit.  When  a  pulse-modulated  RF  signal  is  fed  into  the  cavity,  the  electric 
field  inside  the  cavity  undergoes  a  build-up  and  a  decay  process.  By 
subtracting  the  transient  incident  voltage  pulse  from  the  transient  cavity 
voltage  signal,  we  have  the  transient  reflected  voltage  signal.  The  shape  of 
the  reflected  voltage  pulse  varies  with  p.  Figure  11  shows  three  different 
reflected  pulses  corresponding  to  <  1,  )9  *  1,  and  >  1.  From  transmission 
line  theory, we  have 

P  -  X/(2  -  X)  (6) 


where  X  =  b/a;  b  is  the  height  of  the  latter  peak  in  the  reflected  pulse  and 
a  is  the  height  of  the  earlier  peak  of  the  reflected  envelope. 


The  transient  re¬ 
flected  voltage  of 
the  cavity  was  mea¬ 
sured  at  the  re¬ 
flected  end  of  the 
HP  778D  directional 
coupler  as  shown  in 
Figure  6 .  Corre¬ 
sponding  to  each 
incident  power  the 
reflected  pulse  was 
recorded  and  its 
value  calculated. 

The  reflection 
method  for  measur¬ 
ing  p  is  appropri¬ 
ate  when  the  two 
peak  heights  are 
comparable .  Based 
on  an  error  analy¬ 
sis  we  know  that 
the  measurement 
error  for  X  is  se¬ 
riously  amplified 
when  X  approaches 
2 .  Therefore ,  the 
external  Q  values 
of  the  main  coupler 
ought  to  be  adjust¬ 
ed  to  make  its  cou¬ 
pling  factor  as 
close  to  1  as  pos¬ 
sible  to  minimize 
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Figure  12 .  Measurement  Results 


any  measurement  error.  The  best  way  is  to  use  a  main  coupler  with  an  adjust¬ 
able  antenna.  Figure  12  shows  the  reflected  pulse  during  a  superconducting 
cavity  experiment. 
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(3)  Measurement  of  Q  by  Decay 


Experimentally,  the  loaded  quality  factor  can  be  measured.  With  the 
corresponding  p,  Qg  can  be  calculated.  may  be  measured  by  frequency  scanning 
or  decay.  When  Q  is  as  high  as  10®~10®,  fg  can  no  longer  be  accurately  measured 
by  scanning.  However,  the  decay  technique  is  suitable  for  measuring  a  high  Q 
value.  When  a  stable  electromagnetic  field  is  established  in  the  cavity,  the 
RF  signal  is  cut  off  to  measure  the  decay  curve  of  the  field.  A  time  constant 
T  can  be  determined  from  the  decay  curve.  The  value  of  Qj  can  be  calculated 
from  the  following: 


Qi  -  IT  X  fg  X  r  (7) 

where  fo  is  the  resonant  frequency. 

The  other  curve  in  Figure  12  is  the  free  decay  signal  of  the  electromagnetic 
field  in  the  cavity. 

3.  Control  of  Experimental  Apparatus  and  Data  Acquisition 

In  terms  of  microwave  measurement  technology,  we  have  successfully  located  the 
resonant  frequency  of  the  high-KJ  cavity  and  kept  the  superconducting  cavity 
operating  in  resonance  with  stability  using  a  PLL  technique.  These  methods  are 
new  to  the  microwave  measurement  field  in  China.  In  the  meantime,  micro¬ 
computer  control  and  automatic  data  processing  were  made  possible  for  the 
system  shown  in  Figure  6.  Programs  were  written  for  an  HP  PC-312  microcomputer 
to  control  the  entire  measurement  system  (signal  source,  oscilloscope, 
frequency  counter,  and  power  meter)  to  make  measurement  accurate,  intuitive, 
fast,  and  convenient. 

The  HP  PC-312  has  an  HP  BASIC  controller  which  has  its  own  CPU  on  board.  The 
controller  board  has  a  2.5  MB  BASIC  dedicated  memory  for  remote  control  and 
data  acquisition  through  an  HP-IB  interface.  The  program  written  is  specifi¬ 
cally  for  the  Q  values  of  the  RF  system  and  its  capabilities  include  on-line 
calculation,  screen  (result)  printing,  data  storage,  data  display,  and  post 
data  processing.  A  menu  is  used  to  select  and  modify  parameters.  It  is  simple 
to  operate.  The  software  is  comprised  of  top-down  modules.  It  has  a  solid 
structure  and  has  a  great  deal  of  fault-tolerance.  Its  user  interface  is  also 
friendly.  The  program  can  perform  the  following: 

(1)  Measure  Q  by  way  of  frequency  scanning  or  decay  and  display  results  on  the 
oscilloscope  screen. 

(2)  Print  any  position  of  the  oscilloscope  screen  with  a  laser  printer  in  four 
different  sizes.  Several  figures  may  be  printed  on  the  same  page. 

(3)  Oscilloscope  waveform  data  is  stored  on  disk  for  future  recall,  recalcula¬ 
tion,  and  post  treatment. 

(4)  Redisplay  and  computation  of  stored  waveform. 
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Figure  13 .  The  Function  of  Program  Control 


Figure  13  shows  a 
brief  outline  of  the 
function  of  the  con¬ 
trol  program.  Figure 
14  shows  the  main 
menu  of  the  program. 
Figure  15  is  a  pic¬ 
ture  of  the  1.5  GHz 
Nb  superconducting 
cavity.  Figure  16 
shows  the  assembled 
Nb  cavity.  Figures 
17  and  18  are  photo¬ 
graphs  of  the  micro- 
wave  PLL  measurement 
system  and  the  cryo¬ 
stat  ,  respectively . 
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The  upper  part  of 
Figure  12  shows  the 
calculated  Q  and 
Qo  can  be  derived 
based  on  these  num¬ 
bers.  Program  controlled  automatic  measurement  is  very  convenient  and  timely. 
It  improved  our  efficiency  and  success  rate. 


User  1 

Figure  14.  Main  Menu  of  Program 


VII.  Conclusions 


The  superconducting  cavity  group  at  the  Institute  of  Heavy  Ion  Physics  of 
Beijing  University  has  made  significant  progress  in  the  study  of  the  RF 
superconducting  cavity  in  a  relatively  short  period  of  time.  This  research 
will  serve  as  a  basis  and  technical  foundation  for  future  use  of  the  RF 
superconducting  cavity  and  development  of  superconducting  high-brightness 
electron  injectors. 
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Figure  17.  Microwave  Measurement  System 


Figure  18.  Photo  of  Cryostat 
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[Text]  Abstract 

The  electron  beam  (e-beam)  from  the  wiggler,  after  passing  through  the  180* 
bend,  enters  the  energy  and  time  spectrum  diagnostic  system.  The  180°  bend  is 
designed  for  achromatic  and  isochronous  transport  of  the  e-beam.  Isochronous 
transport  satisfies  the  diagnosis  of  the  time  resolution  of  micropulses  and 
achromatic  transport  satisfies  the  diagnosis  of  the  time  resolution  of  macro¬ 
pulses  .  The  time-dependent  macropulse  energy  spectrometer  is  comprised  of  four 
matching  lenses,  a  two-slit  device,  a  sweep  deflector  coil,  a  90°  symmetric 
bifocal  analyzing  magnet,  a  fluorescent  target  at  the  focal  plane,  a  video 
camera  system  with  an  image  intensifier,  and  an  on-line  computer-based 
graphics  processing  system.  Its  energy  resolution  is  better  than  1,000  and 
time  resolution  is  higher  than  60  ns. 

[Introduction] 

Figure  1  shows  the  layout  of  the  diagnostics  for  the  time-dependent  energy 
spectra  of  micro  and  macropulses  of  the  e-beam  of  the  Beijing  FEL.  Limited  by 
available  space,  the  e-beam  from  the  wiggler  must  be  bent  by  180°  before 
entering  the  time-dependent  energy  spectrometer  system.  This  paper  describes 
the  180°  bend  and  the  macropulse  diagnostic  system. 
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Figure  1.  Beamline  Configuration  for  E-Beam  Diagnostics  of  Beijing  FEL 


I.  Design  of  the  Isochronous  and  Achromatic  Bend 

The  180®  bend  may  be  designed  for  either  isochronous  or  achromatic  transport, 
or  both.  Isochronous  transport  satisfies  the  diagnostics  of  the  micropulse 
time  spectrum  I(t)  and  achromatic  transport  satisfies  the  diagnostics  of  the 
time-dependent  macropulse  energy  spectrum  I(E,t). 

For  a  system  comprised  of  a  horizontal  bending  magnet,  a  quadrupole  lens  and 
a  drift  segment,  its  first-order  six-dimensional  transport  matrix  is  as 
follows^: 
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Furthermore,  when  it  satisfies  the  following  conditions: 


it  is  achromatic.  This  means  that  the  emissivity  of  the  e-beam  does  not 
increase  as  it  passes  through  the  system  despite  the  fact  that  its  energy  is 
dispersed.  When 


fWsi  =  0 

j  W52  =  0  (3) 
=  0 

the  system  is  an  isochronous  transport  system.  This  means  that  the  deviations 
in  trajectory  and  energy  of  the  e-beam  are  independent  of  the  initial  motion 
on  the  bending  plane  and  the  longitudinal  length  of  the  electron  cluster 
remains  unchanged. 

In  engineering  terms,  equations  (2)  and  (3)  imply  that  the  system  must  have  at 
least  five  optical  elements,  such  as  quadrupole  lenses,  to  provide  complicated 
adjustments  and  diagnostics  to  put  the  system  in  isochronous  and  achromatic 
operation.  Obviously,  this  is  very  inconvenient.  To  simplify  the  design  and 
for  ease  of  use,  the  180°  bend  was  designed  symmetrically.  Equations  (2)  and 
(3)  can  be  simplified  and  the  transport  matrix  elements  at  the  center  of  the 
system  (on  the  symmetry  plane)  can  be  made  to  satisfy  the  following  conditions 
(proof  for  isochronous  and  achromatic  transport  is  omitted) . 


=  0  (^) 
"22W51  -  ^21^52  =  0 

_  _  _  _ _  (5) 

In  practice,  two  optical  elements  for  adjusting  the  beam  parameters  can  be 
eliminated.  Considering  the  fact  that  the  streak  camera  used  to  measure  the 
Cerenkov  radiation  of  the  micropulse  is  a  block-type  measurement,  i.e., 
micropulse  and  macropulse  diagnostics  cannot  be  done  simultaneously,  the  180° 
bend  was  designed  to  have  two  pairs  of  S3rmmetric  quadrapole  lenses  (each  pair 
in  series).  When  performing  macropulse  diagnostics,  only  one  pair  is  required. 

This  kind  of  symmetric  bend  has  several  advantages .  For  example ,  when  the 
system  operates  in  the  achromatic  mode,  equation  (4)  becomes; 

[^22£>21  +  +  -^2021)] 

[(^12  ^2^22^  ^21  ^  ^22^22^^26  ^  ^26 

where  H,  jj,  and  Q  are  the  transport  matrixes  of  the  first  bending  magnet, 
second  bending  magnet,  and  quadrapole  lens,  respectively;  Ij  and  I2  are  the 
distance  from  the  exit  of  the  first  bending  magnet  to  the  inlet  of  the  lens 
and  the  distance  from  the  exit  of  the  lens  to  the  inlet  of  the  second  bending 


(6) 

=  0 
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magnet,  respectively.  By  substituting  the  parameters  of  or  system  into 
equation  (6) ,  we  have 

24.787  ksin  kl  -  0.03965  -2.6895cos  kl  -  0.707  -  0 

k 

where  k  =  (G/p)^^^,  G  is  the  gradient  of  the  lens  (Gs/cm)  ,  p  is  the  electron 
magnetic  rigidity  (Gs-cm) ,  and  1  is  the  effective  length  of  the  lens  (1-10 
cm) .  A  unique  solution  of  k  may  be  obtained  from  this  transcendental  equa¬ 
tion.  In  practice,  the  excitation  current  of  the  lens,  Iq,  is  set  based  on  the 
magnitude  of  the  excitation  current  of  the  bending  magnet  (i.e. ,  electron 
magnetic  rigidity),  Ig,  to  make  k  =  kg  which  puts  the  system  in  the  achromatic 
mode.  Furthermore,  the  Iq-Ib  curve  also  serves  as  an  operational  curve  which 
is  precalculated.  When  the  system  is  set  for  the  isochronous  mode,  there  is 
another  set  of  Iq^,  Iq2~Ib  curves.  Similarly,  these  curves  can  be  calculated  in 
advance . 

In  mid-1990,  the  Beijing  FEL  group  requested  that  the  micropulse  measurement 
system  be  installed  behind  the  first  90®  bend.  The  original  achromatic  design 
for  macropulse  diagnostics  was  not  affected  by  this  decision.  However,  the 
symmetric  isochronous  transport  design  would  not  meet  the  requirements  for 
micropulse  diagnostics.  In  operation,  careful  adjustments  of  the  parameters  of 
the  two  lenses  in  the  first  90®  bend  are  required  to  achieve  isochronous 
transport  in  the  90®  bend.  Computation  done  with  the  TEIANSPORT  program^ 
indicates  that  isochronous  transport  is  still  possible  in  this  case  to  meet 
the  requirement  for  micropulse  diagnostics.  Furthermore,  the  calculated 
results  can  provide  guidance  for  operation.  Figure  2  shows  the  calculated 
e— beam  envelope  in  achromatic  and  isochronous  transport.  The  initial  state  of 
the  e-beam  is  that  it  is  in  the  center  of  the  wiggler.  The  e-beam  has  two 
waists  and  the  diameter  of  the  waist  is  1  mm.  The  normalized  emissivity  is 
30?r  mm-mrad.  The  focusing  in  the  y  direction  provided  by  the  B,  component  of 
the  wiggler  field  is  linearized  as  follows : 


y 

y 


cos(z/p,,)  Pj,sin(z/Pp) 
(-1/Pj,)sin(r/Pp)  cos{z/Pp) 


y 

y' 


lz-0 


where  =  N  Ag/2jr. 

N  is  the  number  of  periods  of  the  wiggler  field  and  Ay  =  3  cm. 
II.  Macropulse  Time -Dependent  Energy  Spectrometer 


(7) 


The  system  consists  of  four  matching  quadrupoles  (after  a  180®  bend)  ,  a  double 
slit  device,  a  deflector  sweep  coil,  a  90®  symmetric  bifocal  analyzing  magnet, 
a  fluorescent  target  at  the  focal  plane,  a  video  camera  with  image  intens- 
ifier,  and  an  on-line  computer  graphics  system  (Figure  1).  The  following  is  a 
description  of  the  design  of  major  components. 
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1.  90°  Bifocal  Analyzing  Magnet 


After  achromatic  bending  of  180°,  the  e-beam  is  transformed  into  a  double 
iwalst  (object)  at  the  horizontal  slit  of  the  double-slit  device.  An  image  is 
formed  on  the  fluorescent  target  after  passing  through  the  90°  bifocal 
analyzing  magnet.  This  analyzing  magnet  is  designed  to  be  symmetrically 
bifocal  for  more  stable  operation.  Its  bending  radius  R  =  400  mm,  horizontal 
dispersion  Dh  -  4,  and  amplification  factor  M  -  1.  With  a  full  object  slit 
width  of  1.6  mm,  the  energy  resolution  is  1,000.  This  means  that  10  measure¬ 
ment  points  can  be  obtained  on  the  focal  plane  for  an  e-beam  with  a  ±0.5% 
energy  dispersion.  In  this  case,  the  full  image  width  is  16  mm.  The  designed 
maximtun  e-beam  energy  dispersion  is  ±3%,  which  corresponds  to  a  full  image 
width  of  96  mm. 
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Figure  2.  Envelope  of  Electron  Beam 
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In  order  to  ensure  magnetic  field  uniformity  in  the  working  area  (better  than 
10~^)  ,  the  width  of  the  pole  is  four  times  the  air  gap  between  poles.  The 
operating  magnetic  flux  is  0,l--0.35  T.  In  order  to  ensure  the  same  effective 
hard  magnetic  boundary  at  different  field  strengths,  the  boundary  is  designed 
to  be  a  Rogowski  approximation  boundary  with  magnetic  shield  clamping.  The 
incident  and  emergent  angles  of  the  analyzing  magnet  are  designed  to  be 
tunable  over  the  range  of  25®  ±  15°  so  that  the  actual  incident  and  emergent 
angles  may  be  adjusted  according  to  the  measured  edge  field  distribution. 
Furthermore,  this  arrangement  will  make  it  more  convenient  in  the  future  to 
compensate  for  aberrations  near  the  poles  with  a  circular  arc  boundary 
provided  by  a  hexapole  field.  Figure  3  shows  the  uniform  field  distribution 
calculated  using  the  POISSON  program.^  The  data  (Table  1)  indicates  that  the 
uniform  region  is  over  4  cm  wide. 
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Computation  done  with  the 
TRANSPORT  program  shows 
that  when  the  normalized 
emissivity  of  the  e~beam  is 
SOtt  mm-mrad,  the  object 
slit  can  vary  between  1.2 
and  1.6  mm  without  affect¬ 
ing  the  transport  charac¬ 
teristics  of  the  system. 
The  corresponding  resolu¬ 
tion  of  the  energy  spectro¬ 
meter  is  1,300-1,000. 


Figure  3.  Calculated  Results  of  the 
POISSON  Program 


Table  1. 


X  (mm) 

-30 

-15 

+15 

+30 

AB/Bo 

5.4  X  10"* 

1.6  X  10"'' 

5.7  X  10-5 

5.8  X  10-* 

2.  Sweep  Deflector  Coil 

The  diagnostics  of  the  time-dependent 
macropulse  energy  spectrum  I(E,t)  are 
accomplished  by  a  vertically  sweeping  e- 
beam  generated  by  a  sweep  deflector  coil 
within  a  macropulse.  In  this  case,  the 
vertical  beam  size  on  the  planar  fluo¬ 
rescent  target  represents  the  width  of 
the  macropulse.  The  macropulse  width  of 
the  Beijing  FEL  is  2-5  fis.  Hence,  the 
sweeping  magnetic  field  rises  very  rap¬ 
idly.  In  order  to  avoid  a  very  large 
inverse  peak  voltage  between  the  magnet¬ 
ic  pulsed  power  supply  and  the  excita¬ 
tion  coil,  a  low- inductance  coreless 
saddle-shaped  deflector  coil  is  used  to 
generate  the  uniform  bipolar  magnetic 
field.  The  coil,  whose  structure  is 
shown  in  Figure  4,  is  split  into  two 
sjrmmetric  halves  which  are  tightly 
wrapped  around  the  beam  pipe.  The  angu¬ 
lar  distribution  of  the  number  of  wind¬ 
ings  varies  in  a  manner  close  to  a  co¬ 
sine  function.  Since  the  coil  inductance 
increases  with  the  square  of  the  number 


Figure  4.  Scheme  of  Slow-Def lector 
Sweep  Coil 


of  windings,  the  number  of  windings  should  be  as  small  as  possible.  However, 
the  fewer  the  number  of  windings,  the  poorer  the  uniformity  of  the  bipolar 
field  and  the  larger  the  error  of  the  diagnostics  of  the  time-dependent 
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spectrum.  We  were  required  to  use  the  smallest  number  of  windings  provided 
that  the  uniformity  of  the  field  is  less  than  1  percent.  A  variety  of  schemes 
were  compared  by  computation  using  the  MPVP  program^  and  by  actual  measurements 
using  a  simulator.  The  final  design  chosen  has  a  total  of  52  windings  and  its 
angular  distribution  is  shown  in  Figure  4.  The  conductor  is  a  flat  copper  wire 
doubly  wrapped  with  glass  and  has  a  cross  section  of  1.8  x  2.65  mm.  According 
to  actual  measurements/  the  uniformity  of  the  bipolar  field  is  better  than 
1  percent  when  R  <  15  mm.  The  total  inductance  of  the  deflector  coil  is  257 
/iH. 

If  the  uniform  bipolar  magnetic  field  B  generated  by  the  deflector  coil  rises 
linearly  with  time,  then  the  following  relation  exists  when  the  e-beam  is 
deflected  by  a  small  angle  {  (Figure  5): 


+  2  EE^ 
21(L  +  1/2) 


(8) 


where  E  is  the  kinetic  energy  of  the  e-beam  (MeV) ,  Eq  =  0.511  MeV,  L  is  the 
distance  from  the  exit  of  the  deflector  coil  to  the  fluorescent  target  (m)  ,  Y 
is  the  vertical  coordinate  of  the  scanned  e-beam  on  the  fluorescent  target 
(m)  ,  and  1  is  the  effective  length  of  the  deflector  coil.  Actual  measurements 
showed  that  1  =*=  0.2032  m  and  L  =  1.8783  M.  In  order  to  have  the  image  of  a  35 
MeV  e-beam  fall  at  Y  =  25  mm,  the  deflecting  magnetic  field  should  be: 


B  = 


200>/35^  +  2  X  0.511  X  35 
3  X  0.2032 (1.8783  +  0.1016) 


X  0.025  =  146.3 


Figure  5 .  Relation  Between 
Drive  Current  and  Deflecting 
Power  of  the  Sweep  Coil 
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Figure  6.  Magnet  Excitation  Curve  of 
the  Sweep  Coil 


Figure  6  is  the  measured  magnetic  excitation  curve  of  the  deflector  coil.  It 
is  obvious  that  the  excitation  current  required  is  less  than  30  A.  The  design 
specification  of  the  magnetic  excitation  pulsed  power  supply  is  that  the  pulse 
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current  amplitude  is  ±30  A.  The  leading  edge  rises  linearly  and  the  rise  time 
is  2-5  fis,  adjustable  at  different  scales.  This  power  supply  is  being 
developed.  The  maximum  inverse  peak  voltage  of  the  deflector  coil  is: 

L-^  =  257  X  10-«  X  - =  7  .71  (kV) 

dt  2  X  10-« 

which  meets  the  voltage  breakdown  requirement  of  the  amplifier  of  the  power 
supply. 

Because  the  vertical  "object  slit"  is  designed  to  be  1/6  mm  in  full  width,  the 
effective  working  width  of  the  fluorescent  target  in  the  vertical  direction  is 
50  mm.  Hence,  within  a  macropulse,  it  is  possible  to  obtain  30  points  in  the 
time-dependent  energy  spectrum  for  diagnostics.  For  a  macropulse  which  is  2-5 
tis  wide,  the  camera  is  required  to  have  a  resolution  of  better  than  60-150  ns. 


3.  Fluorescent  Target,  Video  Camera,  and  Computer  Graphics  Processing  System 


An  imported  French-made  Cr-activated  AI2O3  fluorescent  target  (AF995R)  with  an 
effective  working  area  of  100  x  50  mm  was  used.  The  TV  video  camera  used  is  a 
Model  SGB-llN,  which  has  350  vertical  lines  (<0.2  mm)  and  550  horizontal  lines 
(<0.2  mm)  and  a  time  resolution  of  approximately  8  ns.  The  camera  is  equipped 
with  a  microchannel-plate  light  intensifier  and  the  amplification  factor  is 
10,000.  Further,  the  applied  bias  voltage  between  microchannels  can  be  used 
for  high-speed  shutter  control.  The  shortest  gate  control  interval  is  several 
dozen  nanoseconds.  Therefore,  it  can  be  used  to  take  the  energy  spectra  of  the 
e-beam  at  different  intervals  within  the  same  macropulse  in  order  to  study 
imortant  mechanisms  such  as  the  buildup  process  of  the  FEL  oscillation  and  the 
emission  efficiency  of  a  micropulse.  Figure  7  shows  a  block  diagram  of  the 
image  processing  system. 


The  CA5000-P 
video  image 
acquisition 
card  can  be 
directly 
plugged  into 
an  IBM-PC. 
The  follow¬ 
ing  figures, 
tables ,  and 


Figure  7.  Block  Diagram  of  Image  Processing 


data  can  be  provided  on-line  by  the  energy  spectrometer  by  using  the  image 
processing  software  and  FEL  research  software. 


a.  Macropulse  time-dependent  energy  spectra  (current  vs.  energy,  current 
vs.  time,  energy  dispersion  vs.  time,  instantaneous  current  vs. time, 
etc . ) . 


b.  FEL  emission  efficiency. 
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[Text]  Abstract 

An  analytical  study  is  done  on  mode  locking  in  an  PEL  (free  electron  laser) 
and  the  mode -locking  conditions  are  obtained.  It  is  derived  that  due  to  mode 
coupling  the  number  of  active  modes  in  the  cavity  is  reduced  by  a  factor  of  m,., 
the  low  gain  factor  is  increased  by  a  factor  of  me,  and  the  exponential  growth 
factor  goes  up  by  a  factor  of  nic  is  the  ratio  of  the  micropulse 

repetition  frequency  to  the  frequency  interval  of  modes  in  the  cavity  Lv^, 
i.e.  ,mc  =  LUf/LVf.. 

I.  Introduction 

In  an  RF-Linac-driven  IR  short-pulse  FEL,  when  the  micropulse  length  Lj,  is 
comparable  to  the  total  slip  between  the  wiggler  and  the  light  beam  A,  i.e., 
I4,  A,  the  propagation  of  the  e— beam,  i.e.,  the  nonsteady— state  term,  must  be 
taken  into  account  in  the  equation.  The  mode  locking  under  such  conditions  has 
been  investigated  by  numerical  modeling  in  reference  [1] .  However,  when  the 
length  of  the  micropulse,  1^,,  is  still  much  larger  than  the  total  slip  between 
the  wiggler  and  the  light  beam  A,  i.e..  Lb  »  A  (such  as  t^~  10  ps,  Ag-l  /im, 
Lb  =  ^z'^b  ~  0-3  cm,  A  =  NyAg  ~  100  x  10~*cm  =  0.01  cm),  the  slip  effect  between 
the  e— beam  and  the  light  beam  can  be  neglected.  It  can  be  treated  as  a  steady 
state.  The  mode-locking  problem  in  this  case  can  be  studied  analytically.  This 
is  an  attempt  to  analytically  study  the  mode-locking  conditions  and  mode¬ 
coupling  gain  of  an  IR  short-pulse  FEL. 
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II.  Gain  of  Coupled  Mode  IR  Short-Pulse  PEL 
1.  E-Beam  Pulse  Profile 


Let  us  assume  that  the  e— beam  pulse  follows  a  Guasslan  profile,  i.e., 

n(to)  -  nof(to)  (1) 


fu.)  .  ^ 

v'SlT 


(2) 


where  ng,  and  T|j  are  the  density,  longitudinal  velocity,  mean  square  root 

of  the  spatial  length  of  the  pulse,  and  time  width  of  the  pulse,  respectively. 

2.  Optical  Mode  Evolution  Equation 

If  the  e-beam  pulse  length  Is  still  greater  than  the  overall  slip  between  the 
e-beam  passing  through  the  wiggler  and  the  light  beam,  the  problem  can  still 
be  treated  as  a  steady-state  case.  According  to  reference  [2],  the  equation 
for  the  PEL  longitudinal  mode  which  takes  the  pulse  structure  of  the  e-beam 
into  consideration  is  as  follows: 

With  a  linearly  polarized  wiggler,  we  have 


4cYo‘ 


5  = 


c 


(3) 

eE„(z) 

(4) 

eB^ 

w 

(5) 

(6) 

m 

(O  -  Ji  (^) 

(7) 

(8) 

2  (l+a3) 

-  - « 

(9) 

_  2k„Yo 

ii2 

(10) 
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where  Ep,  Wp,  and  kp  are  the  complex  optical  field,  frequency,  and  wave  ntimber 
of  mode  n,  respectively;  e,  m,  and  no  are  the  charge,  laass,  and  beam  density 
of  the  electron,  respectively;  and  k„  are  the  amplitude  and  wave  number  of 
the  wiggler  field;  and  Jq  and  are  the  zeroth-  and  first-order  Bessel 
functions  of  the  first  kind. 

In  an  PEL,  due  to  the  structure  of  the  e-beam,  mode  coupling  can  occur  in  the 
linear  region.  Hence, 


F„(z,to)  =  =  1  -  I.„a„dz  (11) 

m  ^  m 

(12) 

yj-'o  Yo  "jj  '  ' 

By  substituting  equations  (11)  and  (12)  into  (3),  we  have 
da^  =  1  ^ 

OZ  4C7Yn<0„  0  0' 


iw. 


Let 


+  1 


,  (2ic,)2  ,  1 


i«„to 


4cyo(»„p2 


*  2 

X  f^dz'  dz" a  (z") 

JO  JO 

=  -f/.i  to) 

^  2 

^2  =  ^/.2 


(13) 

(14) 

(15) 


By  substituting  equation  (2)  into  (14)  and  (15),  we  have 


-»',oV2V», 


(16) 


-  g- (“»-•» »>’«V2v’, 


(17) 


For  an  IR  short-pulse  PEL, 


-  lO^d/cm) 
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lO'^cm 


(l>  0 

(  "  ^)2  -  10« 
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(18) 
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By  substituting  equations  (17)  and  (18)  into  (13),  we  have 


da„(z) 

dz 


4cYoO)„h* 


m 


—  r<iz'  /'*'  dz"e^‘*'^a„(z") 

®  0)  „  Jo  Jo 


(19) 


Equation  (19)  is  the  mode-evolution  equation  of  the  IR  short-pulse  FEL.  As  we 
can  see,  if  the  following  inequality  is  valid: 


2«2 


<  1 


(20) 


then  it  is  possible  to  couple  several  modes.  From  equation  (20),  the  pulse 
length  of  the  e-beam  is  dependent  on  the  number  of  coupled  modes.  This  is 
because  the  farther  apart  wm  from  Wp  is,  or  the  larger  the  nvunber  of  coupled 
modes,  then  the  shorter  is  required  to  be  to  keep  equation  (20)  valid. 

3.  Number  of  Coupled  Modes 


According  to  reference  [1],  the  number  of  modes  that  can  be  simultaneously 
activated  in  the  resonant  cavity  is 


Av , 
Av , 


(21) 


where  Ai/g  is  the  gain  bandwidth.  It  is  at  least  determined  by  the  uniform-gain 
bandwidth  and  Ai/g  >  i/s/2Njy,  where  i/g  is  the  center  frequency  of  the  light  beam, 
Nj,  is  the  number  of  periods  of  the  wiggler  field,  and  Ai/g  is  the  frequency 
interval  between  modes  in  the  resonant  cavity;  “c/2L(;,  where  is  the 
length  of  the  resonant  cavity  which  is  at  least  equal  to  and  L,j  is  the 
length  of  the  wiggler. 


•••  nig  > 


2jl 

2N„ 


(22) 


Therefore,  thousands  of  modes  may  be  simultaneously  activated  in  the  cavity. 


In  reality,  the  number  of  activated  modes  decreases  with  growing  coherent 
length.  Nevertheless,  this  kind  of  growth  is  limited  by  the  pulse  length.  The 
shorter  the  pulse  length,  the  larger  the  number  of  activated  modes. 


As  described  in  reference  [1],  an  interferometer  inside  the  cavity  (e.g.,  a 
slightly  oblique  etalon)  might  be  used  to  turn  independent  pulses  into 
coherent  light  pulses.  The  number  of  light  pulses  that  might  be  coupled  is 
determined  by  the  niomber  of  independent  coherent  pulses .  The  number  of 
independent  pulses  is  further  determined  by  the  repetition  frequency  of  the 
e^beam  pulse.  Since  we  adjust  the  length  of  the  etalon  with  respect  to  the 
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repetition  frequency  of  the  e-beam  pulses  so  that  the  pulses  in  the  etalon 
overlap  with  the  new  ones,  the  number  of  modes  to  be  coupled  can  be  determined 
by  the  following  equation: 


(23) 


where  Ai/p  is  the  repetition  frequency  of  the  electron  micropulses .  It  is  given 
in  reference  [1]  that  Lg  ”  6  m,  -  1  GHz  and  mg  -  40. 

From  the  number  of  coherent  pulses,  mg  -  40,  the  maximum  difference  between  the 
frequency  of  the  coupled  mode  and  the  given  laser  frequency  Wp  can  be 
calculated.  By  substituting  this  difference  and  the  given  electron  pulse 
length  of  Oj  «  0.1  cm  into  equation  (20),  it  is  found  that  the  equation  is 
valid.  Hence,  me  -  40  is  the  number  of  modes  that  can  be  coupled. 

4.  Low-Gain  Situation 


a„  =  |a„|e^*"  =  const 
=  |a„|e^*-  =  const 

By  substituting  equations  (24)  and  (25)  into  (19),  we  have 


(24) 

(25) 
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(1)  Mode -Locking  Conditions 


«  1 


(26) 

(27) 

(28) 
(29) 


2V* 

|a»l  -  \^n\ 
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(n  -  1,  2,  • • • ,  mg;  m  -  1,  2,  • • • ,  me) 


(30-1-4) 
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Equation  (30)  indicates  that  when  the  phases  (and  amplitudes)  of  me  short-pulse 
modes  of  almost  the  same  frequencies  are  identical,  they  are  coupled  into  a 
single  optical  mode. 

(2)  Gain  of  Coupled  Mode 

Using  equation  (30)  in  (26),  we  have 


Gc  —  nieGi 

^  (2k J 

-  7^, 

(3)  Power  of  Coupled  Mode 

Pc  =  Po  (1  +  Gc)  «  Po  (1  +  mcGi) 

5 .  High-Gain  Situation 
Let 


-  3^,6 


By  substituting  equations  (35)  and  (36)  into  (19),  we  get 


(31) 

(32) 

(33) 

(34) 

(35) 

(36) 


X  exp[-  (<j» J^a%/2v%]  exp  [i  (Aic„-AicJ  z+i  (♦„-<!>„)  z+(a„-a„)  z]  (37) 

”  '^1  [(A;c„+*„)=*+o"J2  a„^  «„  *  o„ 

X  exp[-  ((o  „)^a%/2\%]  exp  [i  (Ak^-Ak^)  z+i  „-<!>„)  z+  (a„-a„)  z]  (3^) 

(1)  Mode-Locking  Condition 


(n  =  1,  2, 
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(2)  Gain  of  Coupled  Mode 


By  applying  equations  (37)  and  (38)  to  (39),  we  get 

{mJS,)  [(Air„+*„)2-e“„] 


[Aif„+*„)2+a*„]2 

From  equations  (40  and  (41),  we  can  obtain  the  following  solutions: 


AJr, 


3  (m<,Go) 
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(3)  Power  of  Coupled  Mode 


Pc  =  Poexp{>/r iltt^Ga)  i/5l„[1-  ( 


A  if. 


3(incGo)i/^ 


)"]} 


(40) 

(41) 

(42) 

(43) 

(44) 


III.  Conclusions 

The  gain  of  the  coupled  mode  of  an  IR  short-pulse  FEL  was  studied  and  mode¬ 
locking  conditions  (equations  (30)  and  (39))  were  obtained.  Due  to  mode 
coupling,  the  number  of  active  modes  in  the  cavity  is  reduced  by  a  factor  of 
me  and  the  low  gain  is  increased  by  a  factor  of  me  (equation  (31)).  The 
exponential  growth  rate  will  rise  by  a  factor  of  me^/®  (equation  (42));  me  can 
be  determined  by  equations  (20)  and  (23). 

References 

1.  D.  Oepts  and  W.B.  Colson,  "Phase  Locking  in  an  Infrared  Short-Pulse  Free 
Electron  Laser,"  IEEE  JOURNAL  OF  QUANTUM  ELECTRONICS,  Vol  26  No  24, 
1990,  p  723. 

2.  Yu  Ming  [0060  2404],  "General  Description  of  Various  FEL  Longitudinal 
Nodes,  "  QIANG  JIGUANG  YU  LIZl  SHU  [HIGH  POWER  LASER  AND  PARTICLE  BEAMS]  , 
Vol  3  No  2,  1991,  p  127. 


38 


Effects  of  Mirror  Coating  Thickness  Nonunifonnities  on  Polarization,  Optical 
Image  Properties  of  Laser,  Beacon  Light 

92FE0442D  Chengdu  QIANG  JIGUANG  YU  LIZI  SHU  [HIGH  POWER  LASER  AND  PARTICLE 
BEAMS]  in  Chinese  Vol  4  No  1,  Feb  92  pp  59-64 

[Article  by  Xiong  Shengming  [3574  0524  2494]  and  Zhang  Yundong  [1728  0061 
3159]  of  the  CAS  Institute  of  Optics  and  Electronics,  P.O.  Box  350,  Shuangliu, 
Chengdu  610209:  "Effects  of  Mirror  Coating  Thickness  Nonuniformities  on 
Polarization,  Optical  Image  Properties  of  Laser,  Beacon  Light";  MS  received  23 
Jul  91,  revised  20  Sep  91] 

[Text]  Abstract 

The  distribution  of  coating  thickness  across  a  0 . 8-m- diameter  mirror  is 
discussed  under  the  geometric  constraints  of  the  coating  chamber.  The  effect 
of  nonuniformity  of  the  metal-enhanced  high -reflective -index  coating  on  the 
polarization  of  the  light  beam  is  presented.  Phase  distortions  associated  with 
the  1.06-^m  laser  and  the  0.5893-/im  beacon  light  are  calculated.  By  means  of 
design  modification,  phase  changes  of  the  same  sign  are  obtained  to  allow  the 
beacon  light  to  achieve  its  expected  objective. 

I .  Introduction 

The  diameter  of  the  mirror  on  the  emitting  end  of  a  high  power  laser  is  D, 
Consider  the  light  beam  following  a  Gaussian  distribution;  the  angle  of 
divergence  at  the  waist  is  ^  ~  A/D  (where  A  is  the  wavelength  of  the  laser) , 
Furthermore,  considering  other  nonlinear  optical  effects  such  as  the  thermal 
blooming  caused  by  the  transport  of  high  power  laser  through  the  atmosphere 
and  excited  Raman  scattering,  the  mirror  diameter  D  ought  to  be  large.  In 
practice,  the  main  mirror  of  a  high  power  laser  is  usually  5-20  m  in  diameter. 
A  large  diameter  is  difficult  and  costly  to  make.  Moreover,  the  preparation  of 
its  coating,  especially  the  uniformity  of  the  coating  thickness,  is  extremely 
critical  and  difficult.  The  nonuniformity  causes  polarization  aberration  and 
wave-front  distortion  which  affect  the  far-field  diffraction  characteristics,^ 

Another  important  issue  of  a  high  power  laser  system  is  that  there  are  many 
factors  causing  wave-front  errors  in  the  optical  path.  The  energy  of  the  laser 
beam  is  dispersed,  which  needs  to  be  compensated  using  the  adaptive  optical 
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technique.  However,  in  order  to  make  this  kind  of  correction  possible,  a 
suitable  light  source,  usually  called  the  beacon  light,  is  required  as  the 
information  carrier  for  the  wave-front  error. 


A  phase  shift  is  produced  as  a  light  beam  is  reflected  by  a  coated  mirror. 
Different  phase  shifts  are  generated  by  light  of  different  wavelengths. 
However,  due  to  nonuniformity  of  the  coating  thickness,  phase  distortions  are 
generated.  At  different  wavelengths,  the  phase  distortions  are  also  different. 
It  is  even  possible  to  have  phase  distortions  of  opposite  signs.  Hence,  the 
opposite  effect  will  be  produced  when  it  is  compensated  by  the  beacon  light. 
Thus,  the  expected  result  cannot  be  achieved. 

Based  on  existing  conditions,  the  coating  thickness  distribution  on  a  0.8-m- 
diameter  near-spherical  mirror  is  calculated  by  taking  the  geometric 
constraints  of  the  coating  chamber  into  consideration.  Furthermore,  this 
distribution  is  optimized.  The  phase  distortions  produced  by  this  t3rpe  of 
metal-enhanced  high-reflectance  coating  series  at  the  laser  wavelength,  i.e., 
1.06  /xm,  and  the  beacon  light  wavelength,  i.e.,  0.5893  fim,  are  analyzed.  By 
modifying  the  design  of  the  coating,  the  signs  of  phase  distortions  at  these 
two  wavelengths  are  consistent. 


II.  Effect  of  Physical  Constraints  in  Coating  Machine  on  Coating  Uniformity 


There  are  two  tjqjes  of  uniformities  concern¬ 
ing  the  coating  of  a  large-diameter  mirror, 
i.e.,  orientation  uniformity  and  radial 
uniformity.  Orientation  uniformity  can  be 
easily  compensated  by  spinning  the  mirror  at 
a  high  speed  via  an  evaporating  ring  such  as 
the  one  shown  in  Figure  1.  Usually,  for  a 
large-diameter  mirror,  from  a  safety  stand¬ 
point,  a  low  spin  rate  is  desired;  hence,  it 
is  necessary  to  use  a  ring-shaped  vapor 
source  to  achieve  excellent  orientation 
uniformity.  However,  radial  distribution  of 
coating  thickness  is  more  difficult  to  con¬ 
trol.  Radial  distribution  of  coating  thick¬ 
ness  depends  upon  three  parameters^:  1)  the  F 
value  of  the  mirror  and  the  nonspherical 
shape  of  the  surface,  2)  the  distance  between 
the  vapor  source  and  the  substrate  H,  and  3) 
the  source  offset  D  from  the  vapor  source  to 
the  center  of  rotation  of  the  mirror  D  (off¬ 
set)  .  With  a  spherical  rotation  fixture,  the 
coating  thickness  at  an  arbitrary  point  P  on 
the  substrate  is 


Rotating 
primary 
Mirror  mirror 


Figure  1.  Ring  Source  Coating 
Geometry  Using  Series  of  Vapor 
Sources  Arranged  to  Simulate  a 
Ring  Vapor  Source  of  Diameter 
2D  Placed  Underneath  a  Rotat¬ 
ing  Primary  Mirror 
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where 


-  •  [  {h^+D^*{>^)  {hcosA+psinA)  -2pD^3in^] 

Hit  [A2+(£,+  p)2]3/2  .  [7i2+(2>-p)2]3/2 

h  ■=  H  -  R  (1  -  cosA) 


(1) 

(2) 


sinA  =  (3) 

In  the  above  equations,  m  is  the  total  weight  of  the  evaporated  material,  ^  is 
its  density,  R  is  the  radius  of  the  spherical  mirror,  and  p  is  the  radial 
radius  from  point  P  to  the  axis  of  rotation.  The  relative  thickness  is  the 
ratio  of  tgp  to  tgoi  which  is  the  thickness  in  equation  (1)  when  p  =  0. 

Figure  2  shows  the  optimal  condition  for  an 
0,8-m“diameter  F/2.0  mirror  to  achieve  coating 
uniformity.  From  this  figure,  one  can  see  that 
in  order  to  achieve  better  than  0.5  percent  in 
uniformity,  the  distance  between  the  vapor 
source  and  the  substrate  needs  to  be  0.9  m  and 
the  offset  is  0.63  m.  The  coating  machine 
should  have  a  height  of  1.3  m  and  a  diameter  of 
1.5  m.  However,  our  equipment  does  not  meet 
these  requirements.  Due  to  cost  considerations, 
modifications  will  have  to  be  made  for  the 
existing  equipment.  Therefore,  it  is  difficult 
to  achieve  the  ideal  coating  uniformity.  Figure 
3  shows  the  effect  of  offset  D  on  the  coating 
uniformity  at  a  known  height  H.  From  the  fig¬ 
ure,  when  H  is  0.65  m,  the  offset  should  be 
0.48-0.5  m  to  achieve  a  uniformity  of  ±1.5  per¬ 
cent.  Nevertheless,  the  diameter  of  our  exist¬ 
ing  coating  machine  is  1.2  m  and  the  maximum 
offset  one  can  attain  is  0.46  m.  Figure  4  shows 
the  coating  thickness  as  a  function  of  height 
with  a  bias  of  0.45  m.  This  figure  shows  that  a 
±2  percent  uniformity  may  be  obtained  at  a 
height  of  0.56-0.64  m. 

In  the  above  discussion,  shields  are  not  used.  In  all  uniformity  calculations, 
all  vapor  sources  are  small-area  sources  and  obey  the  cosine  distribution.  The 
mirror  is  considered  to  be  spherical. 

III.  Initial  Design  of  Coating  Series 

The  main  laser  wavelength  is  selected  at  near  1  pm,  e.g.,  1.06  pm.  The  beacon 
light  might  either  be  scattered  laser  light  reflected  by  the  target  or  a  laser 
light  of  another  wavelength  on  a  relay  satellite.  We  are  using  the  artificial 
beacon  light  produced  by  the  resonant  scattering  of  sodium  atoms  in  the  alkali 
metal  layer  formed  by  falling  meteorites  80-100  km  high  in  the  atmosphere  in 
response  to  an  0.5893  pm  laser  as  proposed  by  French  astronomer  Labyrie.  It  is 


Figure  2 .  Optimized  Coating 
Thickness  Uniformity  for 
0. 8-m-Diameter  F/2.0 
Concave  Mirror 
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Figure  3 .  Coating  Thickness  Nonuni— 
formity  for  0.8-m-Diameter  F/2.0 
Concave  Mirror  for  Various  Offset 
Distance  and  for  a  Fixed  Source- 
to-Substrate  Distance  of  0.65  m 
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Figure  4.  Coating  Thickness  Nonuni¬ 
formity  for  0.8-m-Diameter  F/2.0 
Concave  Mirror  for  Various  Height 
and  a  Fixed  Offset  Distance  of 
0.45  m 


assumed  that  a  reflectance  of  greater  than  99  percent  is  required  at  1.06  urn 
and  that  at  0.5893  fim  it  needs  to  be  higher  than  95  percent. 

The  initial  design  of  the  coating  is 

G/Cr/Ag/Al203/(LH)VAir 

where  L  and  H  are  the  one-fourth-wavelength  thicknesses  of  SiOz  and  TaoO,  at 
the  laser  wavelength,  respectively.  Cr  and  AI2O3  are  bonding  layers.  G  is  the 
substrate  and  the  geometric  thickness  of  the  Ag  layer  is  70  nm. 

The  design  with  a  full  or  half  one-fourth-wavelength  dielectric  film  on  a 
metal  layer  is  selected  because  it  is  easy  to  monitor.  Initially,  the  (LH)^ 
structure  was  chosen  because  it  has  wide  bands  of  reflectance  near  the  two 
wavelengths  of  interest.  Ag  was  selected  because  it  has  a  high  reflectance 
band  ranging  from  visible  to  infrared.  Based  on  this  design,  the  maximum 
incident  angle  of  the  main  mirror  is  15".  The  following  is  a  discussion  of  the 
effect  of  coating  nonuniformity  on  the  polarization  characteristics  (elliptic- 
ity  and  phase  distortion)  at  the  two  wavelengths  of  interest. 

IV.  Effect  of  Thickness  Nonuniformity  on  Ellipticity  of  the  Initial  Coatine 
Design  ® 

When  a  linearly  polarized  light  shines  on  the  main  mirror  at  an  oblique  angle, 
the  reflected  light  is  elliptically  polarized.  In  other  words,  there  is  a 
phase  difference  between  the  vector  component  P  which  is  parallel  to  the 
incident  plane  and  the  vector  S  which  is  perpendicular  to  the  incident  plane. 
Furthermore,  the  perpendicular  components  have  different  amplitudes.  Figure  5 
shows  the  changes  in  reflectivity  and  phase  as  a  function  of  thickness 
nonuniformity  at  the  main  laser  wavelength  and  beacon  light  wavelength.  At 
1.06  im,  the  phase  difference  Ap-As  <  1*.  Even  with  a  4  percent  thickness 
nonuniformity,  the  difference  in  reflectivity  is  less  than  0.02  percent. 
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At  0.5893  /im,  the  phase  difference  is  less 
than  3.5®  and  the  reflectivity  difference 
is  less  than  0.2  percent.  In  conclusion,  at 
an  incident  angle  of  15®,  the  reflected 
light  is  elliptically  polarized  in  a  very 
flat  manner.  A  4-percent  thickness  nonuni¬ 
formity  has  very  little  impact  on  its 
ellipticity. 

V.  Effect  of  Coating  Nonuniformity  on  Phase 
Distortion  of  Reflected  Waves  at  Different 
Wavelengths 

The  characteristics  of  nonuniformity  of 
thin  coating  include  different  optical 
thicknesses  and  phase  shifts  (phase  distor¬ 
tions).  Phase  shift  will  cause  image  dis¬ 
tortion.  In  particular,  it  significantly 
affects  the  far-field  diffraction  charac¬ 
teristics  of  the  system.  Figure  6  shows  the 
phase  distortions  of  the  reflected  waves  at 
the  two  wavelengths  of  interest  as  a  func¬ 
tion  of  coating  thickness  nonuniformity  for 
the  initial  design.  As  one  can  see,  at  1.06 
fim,  a  +4  percent  nonuniformity  only  results 
in  an  -11®  phase  shift  or  a  wave-front 
distortion  of  -A/32.7.  A  1-percent  nonuni- 


0.003 


Relative  nonuniformity  (Z) 

Figure  5.  Maximum  Difference  in 
p  and  s  Components  of  Reflectiv¬ 
ity  and  Phase  for  Initial 
Selected  Design  as  a  Function  of 
Thickness  Nonuniformity  for  15® 
Incidence  Angle 


formity  causes  a  phase  shift  of  —2.4® 
or  a  wave-front  distortion  of  A/150.  At 
0.5893  /im,  the  phase  distortion  is 
larger.  A  +4-percent  nonuniformity 
results  in  a  phase  shift  of  140.28®  or 
a  wave-front  distortion  of  A/2.57.  A  1- 
percent  nonuniformity  causes  a  phase 
shift  of  -8.5®  or  a  wave-front  distor¬ 
tion  of  A/42.3. 

Figure  7  shows  the  radial  distribution 
of  phase  distortion  at  the  wavelengths 
of  interest  for  an  0. 8-m-diameter , 
F/2.0  mirror.  The  calculation  is  based 
on  an  assumption  that  the  large  mirror 
surface  is  close  to  spherical.  The 
distance  between  the  vapor  source  and 
the  substrate  is  0.55  m  and  the  offset 
is  0.45  m. 


s 


Figure  6.  Phase  Distortion  in  Degrees 
for  the  Initial  Selected  Design  at 
the  Two  Wavelengths  as  a  Function  of 
Thickness  Nonuniformity  for  0.8-m- 


Opposite  phase  distortion  signs  cause  Distance  F/2,0  Mirror 


other  problems  as  well.  An  optically 


erroneous  image  at  the  wavelength  may  be  partially  corrected  by  means  of  self- 
adaptive  optical  compensation.  If  the  wave-front  distortion  of  another 
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wavelength  has  the  opposite  sign,  the 
correctibtt  at  one  wavelength  will 
cause  the  further  deterioration  of  the 
wave  front  at  the  other  wavelength.  In 
order  to  simul-tsneously  correct  the 
wave  fronts  at  two  wavelengths ,  phase 
distortions  at  these  two  wavelengths 
ought  to  have  the  same  sign.  To  the 
extent  possible,  the  phase  distortions 
of  these  two  wavelengths  must  be  con¬ 
sistent.  To  this  end,  a  modified  de¬ 
sign  was  adopted. 

VI.  Modified  Design  and  Its  Effect  on  Figure  7.  Phase  Distortion  for  Initial 
Polarization  and  Phase  Distortion  Selected  Design  Along  the  Radial  Dis¬ 

tance  of  0.8-m-Diameter  F/2.0  Mirror 

To  have  the  same  phase  distortion  for  Source-to-Substrate  Height  of 
signs,  numerous  simple  designs  were  0.55  m  and  Offset  Diameter  of  0.45  m 
investigated.  Based  on  theoretical 

calculation,  the  modified  coating  was  designed  as  follows: 


G/Al^OjAg/Al^Oj  (LH)  *-|/Air 

H  and  L  are  the  one-fourth-wavelength  Ta205 
and  Si02  correction  coatings  for  1.17  fim, 
respectively,  and  AI2O3  is  the  bonding 
layer .  The  effect  of  thickness  nonuniform— 
ity  on  the  polarization  characteristics  of 
this  modified  design  is  shown  in  Figure  8. 

It  was  found  that  the  reflectivity  of  this 
design  varies  significantly  with  thickness 
nonuniformity  at  an  0. 5893-/im  wavelength. 

The  two  designs  have  comparable  reflectiv¬ 
ities  at  1.06  nm.  Compared  to  the  initial 
design,  the  reflectivity  of  the  modified 
design  is  slightly  higher  than  a  +4-per- 
cent  nonuniformity.  At  0.5893  /zm,  the 
reflectivity  of  the  modified  design  is 
lower  than  that  of  the  initial  design. 

Figure  9  shows  the  phase  distortions 
produced  by  the  modified  design  at  two 
wavelengths.  It  shows  that  a  4-percent 
nonuniformity  results  in  a  -23*  phase  Figure  8.  Maximum  Difference  in  p 
shift,  or  a  A/T5.6  wave— front  distortion,  s  Components  of  Reflectivity 

at  1.06  pm.  A  1— percent  nonuniformity  Phase  for  Modified  Design  as 

causes  a  phase  shift  of  -4.5*  or  a  A/80  a  Function  of  Thickness  Nonuni- 
wave-front  distortion.  At  0.5893  pm,  the  formity  for  15*  Incidence  Angle 
distortion  is  larger:  a  4— percent  nonuni— 

formity  produces  a  -55*  phase  shift,  or  a  A/6.5  wave-front  distortion  and  a 
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Figure  9.  Phase  Distortion  in 
Degrees  for  the  Modified  Design 
at  the  Two  Wavelengths  as  a 
Function  of  Thickness  Nonuni¬ 
formity  for  0.8-m-Diameter 
F/2.0  Mirror 

VII .  Conclusions 


0.024  iUdial  disUDOt  (m) 
0.4-0.3-0.2-0.I  I  0.1  0.2  0.3  0.4 


CL  ^ 

Figure  10.  Phase  Distortion  for 
Modified  Design  Along  Radial 
Distance  of  0.8-m-Diameter 
F/2.0  Mirror  for  Source-to- 
Substrate  Height  0.56  m  and 
Offset  Distance  of  0.45  m 

1— percent  nonuniformity  results  in  a 
-14.5®  phase  shift,  or  a  A/24.8  wave- 
front  distortion.  The  wave-front  distor¬ 
tions  corresponding  to  different  coating 
thickness  nonuniformities  have  identical 
signs  at  these  two  wavelengths. 

Figure  10  shows  the  radial  distribution 
of  phase  distortion  for  the  modified 
design  at  the  two  wavelengths  for  the 
0.8-m  F/2.0  mirror. 


The  analysis  described  above  indicates  that  when  using  the  modified  design 
(adding  one-fourth-  or  one-eighth-- wave  length  dielectric  coatings  on  top  of 
metallic  silver) ,  coating  thickness  nonuniformity  could  produce  the  same  phase 
distortion  signs  at  two  wavelengths.  At  an  incident  angle  of  15®,  the  modified 
design  does  not  have  much  impact  on  its  polarization  characteristics.  The 
reflectivity  in  the  wavelength  range  of  interest  does  not  decrease.  This  is  a 
suitable  design  for  an  0. 8-m-diameter  mirror.  For  ease  of  monitoring,  all 
dielectric  coating  was  controlled  by  a  single  wavelength  to  minimize  any 
variation  caused  by  dispersion  due  to  the  use  of  different  wavelengths  for 
controlling  different  layers. 
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[Text]  Abstract 

On  the  basis  of  work  done  on  vacumn  synthetic-aperture  beam  control,  the 
effect  of  phase  control  errors  (including  tilt,  piston,  and  focus  error’s)  on 
coherent  beam  combination  is  analyzed  and  discussed.  A  corresponding  physical 
model  is  presented.  In  addition,  numerical  computation  was  done  and  is  illus- 
trated  in  the  form  of  three-dimensional  plots. 

I .  Introduction 

The  use  of  a  synthetic  aperture  as  a  coherent  high-power  laser  emission  system 
has  attracted  a  great  deal  of  interest  in  recent  years. ^  By  controlling  the 
wave-front  phase  of  the  light  beams  from  each  subaperture  of  the  synthetic 
aperture  system,  it  is  possible  to  achieve  coherent  combination  of  the  light 
beams  in  the  far  field.  Nevertheless,  phase-control  errors  usually  occur  with 
wave-front  phase  control.  There  are  several  factors  affecting  the  control  of 
phase  which  leads  to  a  decline  in  performance  of  the  synthetic  aperture  system 
and  lowers  the  far-field  peak  of  the  combined  beam.  We  have  studied  the  beam- 
control  concept  under  ideal  conditions  (i.e.,  without  various  factors 
affecting  the  precise  phase  control  of  the  beam).  A  physical  model  for  the 
alignment,  transition  and  coherent  and  noncoherent  combination  of  the  beams 
was  presented.  Furthermore,  theoretical  analysis  and  numerical  computation 
were  performed. 2. 3  On  the  basis  of  our  previous  work,  this  paper  discusses 
factors  affecting  the  precision  of  phase  control  (e.g.,  piston  error,  tilt 
error,  and  focus  error)  and  their  effect  on  coherent  beam  combination. 
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A  physical  model  is  presented  and  numerical  computation  is  carried  out; 
results  are  illustrated  in  the  form  of  three-dimensional  figures. 

II.  Theoretical  Analysis  and  Numerical  Computation 

To  facilitate  the  theoretical  analysis,  let  us  assume  that  the  laser  emits  a 
planar  beam. 

Uo(xi,  yi)  “  Bpi  (Xi,  yi)  (1) 


where  B  is  a  constant  which  is  there  for  use  in  the  analysis  to  follow.  Let 
B  =  1  at  this  point;  PiCxj,  y^)  is  the  pupil  of  the  light  emission  system  which 
is  defined  as  follows: 


Pi(Xi,  y^) 


x!*yii{d/2)^ 

lxl^yl)>{d/2)^ 


(2) 


where  d  is  the  diameter  of  the  pupil. 


Let  us  assume  that  the  synthetic  aperture  is  composed  of  N  -  1  x  m  apertures 
and  the  center  coordinate  of  each  aperture  is  Xs(i,j),  ys(i,j),and  i-l,2,»««,l 
and  j=l,2,***,m.  The  coherent  combined  far  field  is^: 

1  m  ix-X^)  /Xz 

{x-Xf)  /Xz 

exp{-2itj  +  Z^yji,j)]] 

where  A(x,y,z)  -  exp  (jkz)  •  exp  [i-^)  (x^+y^)]  /  (jXz)  ,  F{  }  represents  the 

Fourier  transform,  k  =  In/X  where  A  is  the  wavelength,  and  (xt.yt)  is  the 
coordinate  of  the  target.  Let  x^  -  y^  -  0  (i.e.,  all  the  beams  from  the 
subapertures  converge  on  the  axial  axis),  then  we  have 


Uix.y)  -  A{x,y,z)F{pi{ri) 


x/Xz 

y/Xz 


exp  { -2it j  [ ( i ,  j ) 


(3) 


Obviously,  from  equation  (3),  the  peak  light  intensity  of  the  coherently 
combined  beam  is : 


X*f=  |p(0,0)|*  =  |F{pi  (r^) }  |=*  •  (Jxm)^  ^  ‘  (4) 

where  Iih  is  the  peak  intensity  of  each  individual  subaperture  system  at 

(0,0, z). 
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In  the  following  analysis  and  computation,  let  N  =  16 
lO^km. 


d  =  lm,  A=-l/im,  z  — 


Figure  1  shows  the  intensity  distribution  of  the  coherent  combination  of  beams 
by  the  synthetic  aperture.  (The  numbers  shown  are  relative  values  )  Figure  2 
shows  the  energy  distribution  curve  in  the  far  field;  Ro  represents  the  Airy 
disc  radius:  E(AIRY  DISC)/E(TOTAL)  -  0.868299E+00  indicates  that  the  energy  in 
the  first  airy  disc  is  86.8299  percent  of  the  total  energy.  R©  =  0.61  x  A  x 
z/b,  which  is  in  units  of  0.61  x  X  z/b  where  A  is  the  wavelength,  z  is  the 
distance  from  the  mirror  to  the  target,  and  b  is  the  effective  pupil  radius  of 
the  phase— control  array. 


Figure  1.  The  Intensity  Distribution 
of  Coherent  Beam  Combination 


1.  Effect  of  Piston  Error  in  Phase 
Adjustment^ 


f  (AIRY  DI5C)/£(TOTAL)=86.8% 
Figure  2.  The  Regional  Energy 
Distribution  in  the  Far  Field. 
E  =  86.8% 


When  performing  wave  front  phase  adjustment  and  control  of  a  divergent  beam 
from  a  subaperture,  a  phase  error  in  the  form  of  exp  (Jk^jj),  where  is 
independent  of  Xj  and  y^,  is  the  piston  error.  Hence,  the  synthetically 
combined  far-field  distribution  becomes: 

x/  Xz 

Upix.y)  =  A{x,y,z)  Fip^ir^)]] 

y/Xz 

1  »  (5) 

exp[-2iii(-^x^(i, j)  +  -^y^(i,j'))]  •  exp(jir<p_jj) 

Consequently,  the  coherently  synthesized  peak  intensity  is 


1^(0, 0)|= 

1 


0) 

0  P  “ 

o|  g  exp(jk<p,p 

(6) 
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i  N^, 


In  general.  g  expCjircpii) 

Therefore,  I  mp  ^  •  Iim 

This  means  that  the  piston  error  would  cause  the  peak  intensity  to  fall. 


Figure  3.  The  Intensity  Distribution  With  the  Piston 
Errors 

^^11* *  Other  »  *0) 


£(AIRY  DlSC)/£(TOTAL)-86.5  % 

Figure  4.  The  Regional  Energy  Distribution 
in  the  Far  Field  With  the  Piston 
Errors 

^.^-0),  £-*86.5% 


Figures  3  and  4  show  the  far-field  light 
intensity  distribution  and  energy  distri¬ 
bution,  respectively,  with  Vii  “  tPiz  “  ^21 

=  ip22  =  0.5  A  and  ip\i  *  0.  These  figures  show  that  piston  error  will  lead  to 
a  decline  of  the  peak  intensity  of  the  combined  beam.  However,  it  has  a 
relatively  small  effect  on  the  energy  in  the  first  Airy  disc. 


2.  Effect  of  Tilt  Error^ 


A  tilt  error  is  a  phase  error  term  in  the  form  of  exp{  [  jk[ci jx(x-Xg(i ,  j ) )  + 
Cijy(y“ys(i» j))]  }  which  is  caused  by  the  tilt  and  vibration  of  the  subaperture. 
In  this  case,  the  optical  field  of  the  combined  beam  can  be  expressed  as 
follows : 

J  jn  x/  \z 

Uj.{x.y)  =  A{x,y.  z,)y‘'V  FiPiir^j) -expUkCii^ix^^-Xg) +jkCijy{y^-yg)]  ]  \ 

In  x/  %,Z 

=  A (X, y,  z, )  g  ^  F{pi  (ri)  exp 


•exp  [-2itj  ( ] 

=  AU,y,z)gg  y{Pi(r,  )}| 


(x-c^j^z)  fXz 
{x-c^jyZ)  /Xz 


•exp  [  -2n  j  ( ( i ,  j )  +  -^  y.{  ij ) )  ] 


(7) 
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It  is  obvious  that  the  peak  intensity  becomes 


\Ur(0,0)\^ 


Figures  5  and  6  are  the  intensity  and  energy  distribution  when  c,jx  “  Cfjy  «  0.5 
X  10“  rad.  Based  on  these  figures,  we  can  see  that  tilt  error  not  only  causes 
the  far--field  peak  light  intensity  to  fall  but  also  results  in  a  reduction  of 
energy  in  the  first  Airy  disc.  The  larger  the  tilt  error,  the  more  apparent 
this  effect. 


Figure  5.  The  Intensity  Distribution  With  the 

Tilt  Errors.  « 0.5 x  I0'» rad  ) 


3.  Effect  of  Focus  Error 


I  2  3 

£(AIRY  D15C)/£(TOTAL)««5.6*/i 

Figure  6.  The  Regional  Energy  Distribution  in 
the  Far  Field  With  the  Tilt  Errors 
X  10'*rad),£  =  85.6% 


When  focusing  the  beam,  if  it  is  not 

accurately  focused,  a  focus  error  term  such  as  the  following  will  be  intro¬ 
duced^: 

exp  {jke,j[xi-xs(i,j))2  +  (yi-y8(i,j))2]/2) 

Hence,  the  far— field  combined  beam  becomes: 


x/Xz 


^p(x.y)  ^  A(x,y.z)j^T  -expUke^ji  (x,-x,)^*  (y^-y,)^) /2]  ] 

Jn  y/X 


=  A(x,y,z)  F^p^(T^)exp(JkeiJTl/2)} 

y/Xz 

•exp{-2iij  [-^Xg(i ,  j)  J 

Let  €.  .=A,  ./z,  then 

i  </  Ar 

Ufix.y)^  A(x,y,z)J^  I  {>,  (''i  )exp(7  A:  A,.r'/ 2  z) } 

ytiz 

■  exp{  -2Kj[-~x,{  uj )  +  ^  y^i  i,j  )]  } 


50 


(9) 


Hence,  the  peak  intensity  of  the  combined  beam  is 


1  m 


F[p^{r^)exp{jitAijri/Xz) } 


I't  3-1 


^  ^  P-I  J-1 


= 


(10) 


Figure  7  is  the  intensity  distribution  of  the  coherently  combined  beam  at  z  *= 
5  X  lO^m  without  any  focus  error  (Afj  =  0).  Figure  8  shows  the  intensity 
distribution  at  z  =  5  x  lO^m  with  a  focus  error,  i.e.,  =-  A12  =  A21  *=  A22  “ 

A31  A32  ==  A41  =  A42  0.1  and  Ai3  *=  A14  =  A23  *  A24  *  A33  =  A34  =  A43  ^  A44  == 

-0.1;  the  physical  significance  is  that  some  of  the  subapertures  are  focused 
in  front  of  the  common  target  and  others  are  behind  it.  It  is  apparent  that 
the  peak  intensity  drops  considerably  in  this  case.  Moreover,  the  larger  the 
focus  error,  the  more  severe  the  impact  on  the  beam  combination. 


Figure  7.  The  Intensity  Distribution  of  Coherent 
Beam  Combination  Without  Focus  Errors. 

(A„»0},j»5x  10*  m. 


Figure  8.  The  Intensity  Distribution  With  the  Focus 

Error.  {  A „  -  A „  =  A,,  =  A ,  -  A  „  -  4.,  -  <J,. 
”A^«-0.l).  i  =  5x|0‘  m. 


III.  Conclusions 

A  detailed  analysis  has  been  done  on  factors  affecting  the  accurate  control  of 
the  phase  of  beam  combination  by  means  of  a  synthetic  aperture;  furthermore, 
numerical  computation  has  been  done.  The  results  show  that  the  peak  intensity 
of  the  combined  beam  drops  considerably  because  of  failure  of  accurate  phase 
control.  However,  it  should  be  pointed  out  that  the  factors  considered  in  the 
analysis  were  assumed  to  be  static,  instead  of  random  and  dynamic  variables. 
As  for  laser  transport  in  the  atmosphere,  it  is  more  meaningful  to  investigate 
the  effects  of  random  and  dynamic  variables,  such  as  turbulence,  on  the 
accuracy  of  phase  control.  We  plan  to  work  in  this  area  in  order  to  establish 
a  realistic  model  to  study  the  effect  of  turbulence  on  phase— control  accuracy. 
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[Text]  Abstract 

The  emission  characteristics  of  velvet  cathodes  have  been  experimentally 
studied  under  long  (50-100  ns)  and  short  (3  ns)  pulse  conditions  with 
corresponding  cathode  densities  of  400  A/cm^  and  7  kA/cm^,  respectively. 

The  current-voltage  curve  was  measured.  The  cathode  emission  pattern  was 
obtained  using  the  Cerenkov  effect.  Scanning  electron  micrographs  of  the 
velvet  cathode  were  obtained  and  the  experimental  structure  of  the  velvet 
cathode  was  measured.  The  emission  mechanism  of  the  velvet  cathode  is 
discussed. 

I.  Introduction 

An  intense  relativistic  electron  beam  (IREB)  has  a  very  high  current  density 
and  a  very  low  temperature.  It  is  a  prerequisite  for  a  number  of  devices  such 
as  electron-beam-pumped  excimer  lasers  and  FELs.  Different  applications  have 
different  characteristic  requirements  for  the  e-beam.  An  e-beam-pumped  excimer 
laser  requires  a  beam  with  very  good  uniformity  and  the  plasma  between  the 
anode  and  cathode  should  have  a  low  closure  speed.  A  Raman  FEL,  on  the  other 
hand,  requires  a  high  brightness  and  low  emissivity. 

These  requirements  can  be  met  by  using  different  diodes  and  cathode  materials. 
Recently,  it  was  proven  that  velvet  is  an  excellent  cold-cathode  material.^  Its 
major  advantages  are  as  follows:  First,  its  start-up  time  is  very  short, ^  i.e., 
its  impedance  can  drop  from  infinity  to  its  characteristic  value  instanta¬ 
neously.  Second,  the  anode-cathode  plasma  closure  speed  is  very  slow  and 
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plasma  expansion  is  almost  unnoticeable .  Thus,  the  current  profile  of  a  diode 

id-l  square  »ave.  Third  "ao„piref  to  oJSfr 

However,  the  emission  mechanism  of  the  velvet  cathode  is  still  not  well 
understood.  To  this  end,  some  experimental  work  is  done  in  order  to  study  the 
emission  mechanism  involved.  ^ 

II.  Experimental  Apparatus  and  Results 

In  this  experimental  work,  three  different  IREB  accelerators  are  used-  their 
pulse  widths  are  120  ns,  50  ns.  and  3  ns.  The  3-ns  device  is  called  the 
Febetron  706.  They  are  essentially  composed  of  three  parts:  1)  An  energy 
storage  unit,  which  is  a  Marx  generator  for  the  120-ns  and  3-ns  apparatus.  The 
ns  apparatus  uses  a  high-voltage  transformer  as  its  energy  storage  unit. 

2)  A  pulse  formation  line,  which  is  a  water  Blumlein  line  for  the  longer-pulse 
oil-based  Blumlein  line  for  the  Febetron  706.  3)  A  diode 
which  has  a  rectangular  e-beam  cross  section  for  the  longer-pulse  apparatus 

istiL''^f'''th^’^  short-pulse  device.  The  character- 

istics  of  these  IREB  accelerators  are  tabulated  in  Table  1. 


Table  1.  Specifications  of  the  IREB  Accelerator  Systems 


Parameters 


Model  IREB  source  (I) 


Output  voltage  (kV) 


Output  current  (kA) 


Output  impedance 
(fi) 


Pulse  width  (ns) 


The  shape  and  area 
of  cathode 


The  results  are  as  follows. 

1.  Measuring  Cathode  Emission  Pattern  Using  Cerenkov  Effect 

Real-time  measurement  of  the  e-beam  current  was  measured  by  means  of  Cerenkov 
radiation  using  an  optical  fiber  or  a  quartz  plate. ^  This  measurement  is  real¬ 
time  and  IS  not  susceptible  to  interference. 

This  experiment  was  carried  out  using  IREB  accelerator  (II).  The  experimental 
set  up  IS  shown  in  Figure  1.  Figure  2  [photograph  not  reproduced]  depicts  the 
emission  pattern  of  the  velvet  cathode.  It  clearly  shows  that  the  velvet 


Camera 


Figure  1.  Schematic  of  Experimental 
Arrangement 


0.2 


cathode  emission  is  more  uniform.  Its 
emission  pattern  consists  of  a  number 
of  bright  and  dark  spots.  The  sizes 
and  locations  of  these  points  seem  to 
follow  a  certain  pattern. 


The 


-Ho.5 


Velvet  fibers 

Cathode  base 


Figure  3.  Practical  Construction  of 
Velvet  Cathode  (Unit:  mm) 


structure  of  the  velvet  cathode 
was  observed  using  a  microscope  and 
is  illustrated  in  Figure  3.  The  vel¬ 
vet  fibers  are  essentially  standing 
straight,  and  are  approximately  20  /im 

in  diameter.  There  are  a  number  of  gaps  at  the  back  side  of  the  velvet  cathode 
and  the  sizes  and  locations  of  these  gaps  also  follow  a  certain  pattern.  We 
measured  the  sizes  of  all  the  bright  and  dark  spots.  They  are  essentially 
consistent  with  the  sizes  of  the  velvet  cathodes. 


2.  Recording  Time -Integral  Current 
Density  Distribution  Using  a  1-mm- 
Diameter-Aperture  Faraday  Cup 

This  measurement  was  done  on  the 
short-pulse  Febetron  706  system, 
whose  cathode  is  circular.  Figure  4 
shows  its  relative  current  density 
distribution.  From  this  current 
density  distribution,  one  can  see 
some  small  nonuniformity.  The  cur¬ 
rent  density  at  the  center  of  the 
cathode  is  smaller  than  that  at  the 
edge.  The  peaks  on  both  sides  cor¬ 
respond  to  the  9-mm  diameter  which 
serves  as  a  measurement  of  the 
effective  cathode  size. 


Radius 

Figure  4.  Relative  Current  Density 
Profile  Obtained  by  Smaller-Aperture 
Faraday  Cup 
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in  of  the  Velvet  Fibers  After  1.000  Runs.  Shown 

in  Figure  5  [Photograph  not  reproduced] 

From  the  photographs ,  it  is  apparent  that  the  tip  of  the  fiber  has  melted  Its 

Sat^th  smooth  and  has  been  stretched.  Another  unique  feature  is 

that  the  velvet  fibers  become  brittle  after  the  experiment.  iSey  break  off 

condi^orr'"^®  ^  carbonized.  However,  they  have  not  yet  become 


III.  Discussion  of  Results 

^ere  are  two  theories  that  might  be  used  to  interpret  the  emission  mechanism 

cfth?l''u  similar  to  an  explosive  emission 

.  e  believe  that  it  might  explain  the  mechanism  in  our  experiment  with 

below  lOVcm  which  is  the  low-field  case.  Hence,  electrons  are  initially 
emitted  from  the  tip  of  the  fiber.  Since  the  velvet  is  an  Insulator,  the  tip 
becomes  charged  after  emission.  ” 


For  an  explosive  cathode,  once  the 
applied  voltage  exceeds  the  break¬ 
down  voltage,  the  tips  begin  to 
explode.  This  moment  is  a  function 
of  the  applied  voltage  and  the 
cathode  material  and  is  called  the 
breakdown  delay  time.  In  the  early 
stage  of  discovery,  the  breakdown 
delay  times  for  a  number  of  materi¬ 
als  were  quantitatively  measured. 
Russian  scientists  have  done  a 
great  deal  of  work  in  this  area  and 
the  values  are  shown  in  Figure  6.^ 
The  breakdown  delay  time  is  the 
time  after  a  pulse  voltage  is  ap¬ 
plied  to  a  planar  electrode  to  the 


Figure  6.  Experimental  Breakdown  Delay 
Time  as  a  Function  of  Applied 
Field  and  Electrode  Materials 


measured  the  breakdown  delay  times  of 
Figure  6^  applied  voltages  and  the  results  are  also  included  in 


For  most  emission  materials,  the  following  equation^  applies: 

j|  X  -  0.55[{nVi)/Q„]  (1) 

tj  density  of  the  emitter  (before  its  destruction)  in 

nits  of  A/cm  ,  ti  is  the  breakdown  delay  time  in  seconds;  and  p,  X  and  On  are 

Tiletr  conductivity,  and  resistivity  of  4e  mat^ikl  respec! 

Because  of  the  huge  difference  in  resistivity  between  metal  and  insulator,  of 
the  order  of  a  dozen  orders  of  magnitude,  t^  for  velvet  is  very  short.  This  is 
he  reason  why  the  velvet  cathode  has  a  very  short  "closure  time." 
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Another  unique  feature  of  the  velvet  cathode  is  that  its  plasma  moves  at  a 
very  slow  speed.  Theoretically,  the  velocity  of  the  metal-cathode  plasma  can 
be  treated  as  a  gas  kinetics  problem  dealing  with  adiabatic  diffusion  of  gas 
into  vacuum.  The  velocity  can  be  expressed  as  follows: 

V  =  {[47/(7  -  1)]  X  fj  X  (2) 

where  7  =  Cp.Cy  which  is  the  adiabatic  coefficient,  E  is  the  latent  heat  of 
sublimation  for  the  metal,  and  ij  is  the  overheat  coefficient.  Since  the  velvet 
cathode  plasma  temperature  is  lower  than  that  of  a  metal  cathode,  the  velvet 
cathode  plasma  moves  at  a  much  lower  velocity  than  that  of  a  metal  cathode. 

How  are  electrons  replenished  from  the  metal  substrate  to  the  tips  of  the 
velvet  cathode?  There  is  a  surface  internal  network  theory.  First,  there  are 
field-effect  electrons  in  the  high-field  region  near  the  metal-vacuum 
insulator  interface.  Then,  these  primary  electrons  bombard  the  surface  of  the 
insulator  to  create  secondary  electrons.  Under  the  influence  of  the  external 
field,  electrons  hop  along  the  surface  of  the  insulator  to  form  a  network  on 
the  velvet  surface.  Electrons  are  replenished  to  the  tips  by  this  mechanism. 
A  detailed  description  of  this  mechanism  is  available  in  reference  [9]. 

Recently,  R.V.  Latham,  et  al.,  presented  an  FIHEE  (field-induced  hot  electron 
emission)  model®’^’®  to  explain  the  prebreakdown  electron  emission  mechanism  on 
the  high-voltage  vacuum— insulator  electrode  surface.  We  believe  that  it  might 
be  used  to  explain  the  emission  mechanism  for  the  velvet  cathode  with  short 
pulses  (3  ns). 

With  an  applied  field  between  the  anode  and  cathode,  velvet  is  an  excellent 
dielectric  material.  Hence,  the  external  field  could  penetrate  the  velvet  and 
reach  the  metal  substrate  to  cause  the  energy  levels  to  change. 

There  is  a  critical  field  in  an  insulator,  approximately  10®-10®V/cm.  When  the 
applied  field  exceeds  this  strength,  electrons  can  tunnel  from  the  metal 
substrate  to  the  conduction  band  of  the  velvet.  This  heats  the  velvet  up  and 
accelerates  the  ionization  process  which  causes  an  avalanche  in  the  conduction 
band  of  the  velvet. 

Let  us  assume  that  AV  is  the  voltage  drop  across  the  velvet  fiber,  V  is  the 
applied  voltage  between  the  anode  and  cathode,  d  is  the  distance  between  the 
two  electrodes.  Ad  is  the  thickness  of  the  velvet  and  c  is  the  relative 
dielectric  constant  of  the  velvet,  then 

AV'  -  (-i)  (-^)  V  (3) 
€  a 

In  our  experiment,  c  ~  4,  V  =  300--600  kV,  d  =  8  mm,  Ad  =  1,5  mm.  Here,  AV  — 
15-30  kV,  and  E  -  10^-2  x  lO^V/cm.  Under  such  conditions,  dielectric  materials 
such  as  velvet  have  some  conductive  paths  to  convert  the  velvet  from  a 
nonconductive  to  a  conductive  state. 
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In  the  presence  of  a  high  field,  the  Poolo— Frenkel  effect  plays  a  major  role 
(i.e.,  the  field  enhances  thermal  ionization  between  donors  and  traps).  Then, 
the  overall  conductivity  is  a  function  of  the  applied  field.  In  this  case,  j 
and  AV  may  be  expressed  as  follows: 


J  =  Joexp[(e/kD  (eAV/eegd)^/*] 
or  =  Joexpt  (e/2/cD  ( eA  V/eCoCf)  ^^*3 


(4) 


depending  upon  the  location  and  characteristics  of  the  trap. 


In  the  short— pulse  experiment, 
voltage  and  current  profiles  were 
measured  after  30,  80,  and  170 
runs.  From  equation  (3),  we  have 
the  (AV)^/2  versus  Inl  curve  (Figure 
7).  The  results  are  in  good  agree¬ 
ment  with  the  theory  described 
above . 

In  this  case,  the  current  density 
across  the  velvet  is  very  high 
(approximately  7  kA/cm^) ;  there¬ 
fore,  some  fibers  are  burnt.  How¬ 
ever,  the  entire  pulse  lasts  a  very 
short  time  and  burning  might  have 
occurred  at  the  late  stage  of  emis¬ 
sion. 


Figure  7 .  The  Experimental  Character¬ 
istics  of  Inl  Versus  (AV)^^^ 
of  Velvet  Cathode 
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[Text]  Abstract 

The  physical  mechanism  of  microwaves  excited  by  two  relativistic  electron 
beams  (REBs)  in  a  cylindrical  waveguide  is  studied  based  on  magnetic-fluid 
theory.  A  dispersion  relation  has  been  derived  in  the  presence  of  an  infinite 
axial  guiding  magnetic  field.  An  expression  for  the  growth  has  also  been  found 
at  the  Cerenkov  resonant  point, ^  i.e.,  w  =  kv^^g.  The  dependence  of  growth  rate 
upon  related  parameters  is  also  briefly  discussed. 

I.  Introduction 

There  are  several  types  of  high-power  microwave  generators, ^  such  as  the 
dielectric  Cerenkov  microwave  generator,  virtual  cathode  oscillator  (vircator) 
microwave  generator,  and  beam— plasma  interaction  microwave  generator.  Each 
type  of  microwave  generator  has  its  own  unique  advantages  and  limitations  in 
terms  of  band,  power,  and  structure. 

A  novel  microwave  generator,  i.e.,  a  dual-REB  microwave  generator,  is 
presented  in  this  paper.  Physically  speaking,  when  two  electron  beams  move 
toward  each  other,  they  produce  a  two-stream  instability . ^  if  the  two  beams  are 
placed  in  a  cylindrical  waveguide,  this  instability  can  be  coupled  to  yield 
waveguide  electromagnetic  modes  to  consequently  generate  a  microwave  output. 
Since  the  parameters  of  the  two  electron  beams  can  be  adjusted,  it  is  easy  to 
achieve  resonance. 
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II.  Basic  Eqtiations 

Assume  the  initial  velocities  of  the  two  cold  REBs  are  v^q  = 

-v^o®z»  respectively.  The  two  beams  are  distinguished  by  subscripts  a  and 
The  continuity  equation  is: 

+  V  •  =0,  (i  =  a,  P)  (1) 

The  conservation  of  momentum  equation  is 

=  -e{E+  —  Vi  X  B)  ,  a  =  a,  P)  (2) 

d  t  c 

where  P,-  -  m7,-Vj  represents  the  momentum  of  the  electron;  m,  v,- ,  and  -e  are  the 
electron  mass,  velocity,  and  charge,  respectively;  c  is  the  speed  of  light  in 
vacuum;  nj  is  the  electron  density,  and  7i  “  (1  -  is  the  relativistic 

factor  of  the  electron.  E  and  B  represent  the  electric  and  magnetic  fields, 
respectively,  and  satisfy  the  Maxwell  wave  equation 

(V*  -  -L  £  »  -ii  +  V(V  •  E)  (3) 

where  ^  (-ev,-ni)  is  the  total  current  density. 

Equations  (1),  (2),  and  (3)  are  the  basic  equations  of  the  magnetic -fluid 

model . 

III.  Dispersion  Relation 

When  an  infinitely  large  axial  guiding  magnetic  field  is  applied  to  a 
cylindrical  waveguide,  the  two  REBs  only  move  in  one— dimens ion  along  Let 
us  assume  that  the  two  REBs  have  the  same  radius  r^,  and  the  radius  of  the 
waveguide  is  R,  the  mean  densities  of  the  two  beams  are  n^g  I'SO  •  respec¬ 
tively.  The  dispersion  relation  may  be  derived  from  perturbation.  To  this  end, 
we  assume  that 


and 


Vi  -  Vio  +  Vfi,  ni  -  nio  +  nil,  E  “  Ei.  5  ”  Bi 


(Vii,  nil,  El. 


Bi)  oc 


(4) 

(5) 


where  the  subscript  "1"  represents  a  first-order  perturbation  term,  and  w  and 
k  are  the  angular  frequency  of  the  electromagnetic  wave  and  wave  number  along 
Cl,  respectively. 


The  two  REBs  in  one-dimensional  motion  inside  a  cylindrical  waveguide  can  only 
interact  with  the  TMgn  mode,  under  the  following  boundary  conditions: 


^,lr. 


0,  ^,1 


r-xj, 


^,lr.x 


^1  -  ^£*1 

dz  ■  dr 


(6) 
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The  following  dispersion  equation  can  be  derived  from  equations  (1)~(3): 


Yq  (pR)  jj  (pr^,)  -  Jq  {pR)  yI  (px^y  _  Ji(k^x^) 
Y^(pR)Jo(pr^)  -  Jo{pR)Yo{pZi,)  ^  p  Joik^Tj,) 

where 

kl  =  [1-  __J?£5 - p2 

Y.o(ifV-,o-<o)2  Ypo(icvpo+»)2 


(7) 

(8) 


(9) 


Jo  and  Yo  represent  the  zeroth-order  Bessel  function  of  the  first  and  second 
kind,  respectively.  In  equation  (8) 


“pi  =  y,o  =  (1  -  i  =  a,  P 


IV.  Growth  Rate  of  Instability 


The  dispersion  equation  (7)  may  be  solved  numerically.  Let  Aon/R  be  the  nth 
root  (n  =  1,2  -  •  •)  of  in  equation  (7)  .  When  R  =  rt,  Aon  is  the  root  of  Jo(Aon) 
=  0.  Hence,  from  equation  (8),  we  have 


Y«o(<^-kv,o)2  Y|o(‘»+kv.(,) 


-]  (u*  -  k^c^)  = 


Equation  (10)  clearly  points  out  that  the  physical  process  in  a  dual-REB 
microwave  generator  is  the  process  that  leads  to  the  coupling  of  the 
electromagnetic  mode  of  the  waveguide  due  to  two-stream  instability. 


If  both  beams  have  similar  electron  velocities  and  densities,  it  is  possible 
to  find  the  growth  rate  of  the  pure  growth  mode  where  the  real  part  of  the 
frequency  is  close  to  zero  from  equation  (10) .  Since  the  frequency  is  near 
zero,  this  mode  is  of  no  interest  to  us.  Hence,  let  us  consider  the  following 
situation  where  the  density  of  electron  beam  a  is  much  smaller  than  that  of 
electron  beam  y9,  i.e.. 


<0^ 

c 

Y«o 

Ypo 

(11) 


Near  the  Cerenkov  resonant  point  w  -  kv„o,  because  there  might  be  some  strong 
interaction  between  the  electromagnetic  wave  and  electron  beam  a,  the  maximum 
growth  rate  of  instability  can  be  obtained  from  equation  (10) . 
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1/3 


r  = 


2*'^y,o 


[(ySo-1)  (1+^)+1] 

^mo 


(On 


(12) 


where 


(On  = 


<^pP 


^ojiYao'iio 


y|o(1-^)* 


'•0 


(13) 


The  frequency  of  the  microwave  is 

0)  =  Wo  -  T/y/T  «  Wo 


(14) 


V.  Conclusions 

A  magnetic-fluid  model  is  used  to  study  the  interaction  between  two  electron 
beams  and  a  waveguide.  Under  the  approximation  that  the  densities  of  the  two 
electron  beams  differ  significantly,  an  expression,  equation  (12),  for  the 
maximum  growth  rate  near  the  Cerenkov  resonant  point,  w  **  kVjKO,  is  obtained. 
The  results  indicate  that: 

(1)  The  output  microwave  frequency  is  a  function  of  v„o,  v^o>  AkO.  ^ 

rb.  By  properly  adjusting  some  parameters  while  keeping  others  unchanged,  it 
is  possible  to  achieve  some  degree  of  resonance.  It  might  be  more  effective  to 
adjust  the  electron  density  n^o  of  REB  and  the  electron  velocity  Vo,o  of  REB 
a . 

(2)  The  growth  rate  can  be  increased  by  raising  the  n^o  and  n^o  of  REBs. 

(3)  The  growth  rate  can  be  lowered  by  increasing  the  velocity,  or  energy,  of 
the  REBs. 


(4)  With 


“pp 


3  /I  .  ^^1)2 


Ypo(l  + 


when  n  increases ,  Aon  goes  up  and  the  growth 


rate  of  instability  falls.  Therefore,  the  growth  rate  of  the  lowest-order  mode 
TMqi  should  be  the  highest. 


(5)  In  order  to  illustrate  that  the  electromagnetic  radiation  emitted  by  the 
dual— REB  microwave  generator  is  in  the  microwave  range ,  a  numerical  example  is 
given  below. 


Assume  REB  a  and  REB  P  have  the  same  velocity,  i.e.  ,  70.0-7^0“  1.196  and  v^q 
-  Vflo  “  0.548c.  However,  the  densities  are  quite  different:  n^o  “  2.688  x 
10“cm"3  and  n^o  “  ^-30  x  IQi^cm-^.  The  radius  of  the  cylindrical  waveguide 
R  =  1.5  cm.  The  electron  beams  almost  fill  the  waveguide  (rj,"'R).  From  eq.  (13) 
we  can  calculate  that  Wq  “  3.178  x  10^°s  From  Wq  ”  (27r/A)Vj|,o ,  the  wavelength 
of  the  electromagnetic  wave  was  calculated  to  be  3.248  cm. 
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